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Abstract 

Non-relativistic physics is often associated with atomic physics and low-energy phenomena of the strong interactions 
between nuclei and quarks. In this review we cover three topics in contemporary high-energy physics at or close 
to the TeV scale, where non-relativistic dynamics plays an important if not defining role. We first discuss in detail 
the third-order corrections to top-quark pair production in electron-positron collisions in the threshold region, which 
plays a major role at a future high-energy e^e~ collider. Threshold effects are also relevant in the production of heavy 
particles in hadronic collisions, where in addition to the Coulomb force soft gluon radiation contributes to enhanced 
quantum corrections. We review the joint resummation of non-relativistic and soft gluon effects for pair production of 
top quarks and supersymmetric particles to next-to-next-to-leading logarithmic accuracy. The third topic deals with 
pair annihilation of dark matter particles within the framework of the Minimal Supersymmetric Standard Model. Here 
the electroweak Yukawa force generated by the exchange of gauge and Higgs bosons can cause large “Sommerfeld” 
enhancements of the annihilation cross section in some parameter regions. 

Keywords: top quark production, perturbative QCD, (P)NRQCD, ILC, LHC, dark matter 


1. Introduction 

Non-relativistic particle physics is often associated 
with atomic physics, the interactions between nuclei, 
and the low-energy phenomena of the strong interac¬ 
tions of the charm and bottom quarks, which form 
quarkonium bound states. But non-relativistic dynam¬ 
ics also governs the interactions of weak-scale particles 
such as the top quark or the hypothetical supersym¬ 
metric partners of the Standard Model (SM) particles, 
when they are slowly moving. In this review we 
cover three topics in contemporary high-energy physics 
where the unique non-relativistic dynamics caused by 
instantaneous, potential interactions plays an important 
if not defining role; the production of top-quark pairs 
in electron-positron collisions in the threshold region, 
which can provide a measurement of the top-quark mass 
with unchallenged precision and could be realized at 
a future high-energy e^e~ collider. The hadronic pair 


production of top quarks at Tevatron and LHC, or yet 
unobserved heavy particles at the LHC, where in addi¬ 
tion to the Coulomb force soft gluon radiation contrib¬ 
utes to enhanced quantum corrections, which should be 
summed. And finally, the pair annihilation of dark mat¬ 
ter particles within the framework of the Minimal Su¬ 
persymmetric Standard Model (MSSM). While the first 
two processes are determined by the long-range colour 
Coulomb force, the low-velocity annihilation of heavy 
neutralinos in the MSSM can be dramatically enhanced 
by the short-range electroweak Yukawa force generated 
by the exchange of gauge and Higgs bosons. 

We discuss the phenomena and results but also em¬ 
phasize theoretical methods and techniques. Since there 
is always a small relative velocity v in the problem, non- 
relativistic systems involve (at least) the three scales m 
(mass), mv (momentum), tmP' (energy), and it is appro¬ 
priate to formulate effective Lagrangians to organize the 
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calculation, even if the underlying theory — QCD, the 
SM, or the MSSM — is known. The three topics dis¬ 
cussed here have further in common that the relevant 
Lagrangian couplings are weak at all three scales. This 
allows for a systematic treatment with approximations 
of well-defined parametric accuracy, even though ordin¬ 
ary perturbation theory in the coupling breaks down, as 
will be discussed. 

Each of the three topics provides its own specific 
challenges. In top-quark pair production near threshold 
in collisions, it is the high precision of the en¬ 
visioned mass measurement, which demands calcula¬ 
tions with unprecedented third-order accuracy in the 
non-relativistic approximation that can be achieved by 
a combination of effective field theory and multi-loop 
techniques. In hadronic production the required pre¬ 
cision is less, but the coloured partonic initial state of 
quarks and gluons implies that soft-gluon resummation 
must be integrated into the non-relativistic expansion, 
and a more general treatment of colour is needed for fi¬ 
nal states which can be produced in various colour rep¬ 
resentations. Dark matter (DM) annihilation involves 
yet another issue. The non-relativistic “Sommerfeld” 
enhancements from electroweak gauge boson exchange 
are operative only, when the DM particles are in the TeV 
range, such that the inverse Bohr radius scale is of or¬ 
der of or larger than the mass of the electroweak gauge 
bosons. For such heavy, weakly interacting particles, 
the mass splitting between the members of the elec¬ 
troweak multiplet are small, so that co-annihilation ef¬ 
fects at dark matter freeze-out are generic. While the 
non-relativistic treatment is only performed at leading- 
order, the complexity of the problem arises from a large 
number of interacting two-particle states, kinematic¬ 
ally closed channels in the Schrodinger problem, which 
must be solved numerically, and the large number of fi¬ 
nal states that contribute to the pair annihilation cross 
section in a model such as the MSSM. 

The outline of this review is as follows: in the next 
section we briefly summarize several methods and tech¬ 
nical tools, which are necessary to perform the calcu¬ 
lations. In particular, we briefly describe the threshold 
expansion, introduce NRQCD and PNRQCD, and men¬ 
tion a few details in connection to the underlying multi¬ 
loop calculations. Section is devoted to the top- 
quark pair production cross section near threshold in 
annihilation. We introduce the main ingredients 
in the third-order calculation and discuss the numer¬ 
ical results. In Section top-quark and supersymmet¬ 
ric pair production is considered at hadron colliders. 
After describing the joint soft and Coulomb resumma¬ 
tion formalism, the top-pair invariant mass distribution 


near threshold is considered. We then summarize res¬ 
ults on the resummed total top and SUSY pair produc¬ 
tion cross sections, and compare it to fixed-order res¬ 
ults, and results with soft-gluon but without Coulomb 
resummation. Section|^leads through the effective held 
theory treatment of Sommerfeld enhancement in dark 
matter annihilation with a short discussion of a method 
to determine the enhancement reliably in situations with 
several channels without extreme degeneracies. We 
also highlight a few results for a wino-like dark mat¬ 
ter particle and a series of models that interpolates from 
a Higgsino- to wino-like neutralino. We conclude our 
review with a summary in Section]^ 

2. Methods and techniques 

Before addressing more specifically the three topics 
of this review, we summarize the main methods and 
techniques. 

2.7. Threshold expansion 

Close to the production threshold of two heavy 
particles (for simplicity, with equal masses) there are 
three characteristic scales: the mass of the particle, m, 
its three-momentum of order mv, and the kinetic energy 
^fs - 2m = mv^. For example, for top quarks we have 
m xi 170 GeV, mv =» 20 GeV and mv^ ^ 2 GeV, and 
thus a strong hierarchy among the scales. Furthermore, 
mv^ » Aqcd, which means that the strong coupling 
Osip) is always in the perturbative regime. In this re¬ 
view, we only consider systems where this condition is 
satisfied. 

The presence of the small parameter v leads to a 
breakdown of the standard perturbation expansion in 
Qfj due to kinematic 1/v enhancements. A reorganized, 
non-relativistic expansion must be used, where both as 
and V are simultaneously considered as small but as/v 
of order one. In terms of Feynman diagrams, a summa¬ 
tion of certain terms to all orders in Os is required. At 
the same time, since v is small, one does not need the 
full dependence on v of every fixed-order diagram. 

For a given Feynman diagram the expansion in v can 
be constructed without first computing the full expres¬ 
sion using the threshold expansion m. The method ex¬ 
ploits that every diagram can be decomposed into a sum 
of terms, for which each loop momentum is in one of 
the following four momentum regions: 


hard (h) : 

f ~m. 

( ~ m 


soft (s) : 

7° ~ mv. 

( ~ mv 


potential (p) : 

7® ~ mv^. 

{ ~ mv 


ultrasoft (us) : 

7" ~ mv^. 

( ~ mv^ 

(1) 
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In each region the integrand of the Feynman diagram 
is expanded in the small parameters of the respective 
region. Afterwards the loop integration over the com¬ 
plete £/-dimensional space-time volume is performed. 
The property of dimensional regularization that scale¬ 
less integrals vanish, guarantees that no double count¬ 
ing occurs. The arguments in favour of this construc¬ 
tion provided in HI rely on assumptions on the location 
of singularities in the Feynman integrand and examples. 
Further justification can be found in ||2l. It may happen 
that the separation of the integrand into regions gener¬ 
ates spurious ultraviolet and infrared divergences. How¬ 
ever, they cancel in the sum of all contributions. 

2.2. Non-relativistic effective theory 

The procedure described in the previous subsection is 
essentially equivalent to the explicit construction of ef¬ 
fective non-relativistic Lagrangians within dimensional 
regularization, since the terms that arise in the threshold 
expansion can be interpreted as modified propagators 
and vertices. In a certain sense, the threshold expansion 
defines these Lagrangians in dimensional regularization 
by providing the precise rules for performing consistent 
matching calculations. 

To be definite, we refer to a system of a heavy quark 
(Q) and antiquark interacting via gluons. The effective 
Lagrangians are constructed in two steps following 
the following scheme for integrating out momentum 
modes: 

-CQCD[Q{h,s,p), g{h,s,p,us)] p>m 


-Cnrqcd [Q(s, p), g{s, p, Mi)] mv <p< m 


-CpNRQCD [Q(p), g(Mi)] jj. <mv 

In square brackets the modes of the heavy quarks (Q) 
and massless particles (g) are displayed which are still 


contained in the effective Lagrangian; the others are in¬ 
tegrated out when the energy scale p is lowered as in¬ 
dicated on the right. 

The first step leads to NRQCD |[3}21. Its Lagrangian 
has the following structure 

-Cnrqcd = + ■Cx + -Ci/^t- -i- -Cg + -Ciight, (2) 

where explicit expressions of the individual contribu¬ 
tions can be found in 0. The bilinear heavy-quark 
Lagrangians X,/, and contain the kinetic terms for the 
two-component quark and antiquark fields if/ and x and 
the interactions with the chromomagnetic field in an ex¬ 
pansion up to including terms of order 1/m^. con¬ 
tains four-fermion quark-antiquark terms and Xg is the 
pure gluon part, again expanded in 1 /m. Finally, Xiight is 
the same as the light-quark Lagrangian in full QCD. The 
Feynman rules derived from the Lagrangian Q can, 
e.g., be found in 0. All matching coefficients relevant 
for a third-order calculation can also be found in 0. 
A subtlety that is explained there in detail is that the 
short-distance coefficient functions of the operators in 
the NRQCD Lagrangian must be kept (f-dimensional in 
dimensional regularization. 

In the second matching step in the above scheme soft 
modes and potential massless modes are integrated out. 
It has been suggested in the context of the threshold 
expansion in El and at the level of an effective Lag¬ 
rangian in Q. The result is the so-called potential 
NRQCD (PNRQCD) Lagrangian developed in various 
forms in CHni. It only contains ultrasoft light fields 
and potential heavy quarks, which simplifies the scal¬ 
ing in the velocity expansion. Note that since the com¬ 
ponents of the hard momentum integrated out in the 
first matching step are larger than those of the modes 
in the effective theory, the NRQCD Lagrangian is local. 
This can no longer be expected for PNRQCD. However, 
since only the three-momentum of the potential modes 
Qip) is of order of the soft scale mv, the non-locality of 
the PNRQCD Lagrangian refers only to space but not to 
time. The relevant terms in the PNRQCD Lagrangian 
are given by 0 


XpNRQCD - tA p'i9o -I- g.5Ao(f, 0) -H 


2m 8m^ 



/5o + gsAoit, 0) - 


^ \ 

2m 


/■ 


ld-\. 




(x -H r) Vab-cdir, d) 


XcXd 


{x) - gs^Hx)^ ■ E(f, 0)i/r(.x) - gsXXx)x ■ E(f, Q)x{x ), 


(3) 
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where 

Vab-cd(r, d) = + 6Vab-cd(r, d) . (4) 

Vo = —cbslr is the tree-level colour Coulomb poten¬ 
tial. The first two terms of -Cpnrqcd contain the kin¬ 
etic terms including the relativistic corrections and the 
third term is responsible for heavy-quark potential inter¬ 
actions. Note that the heavy-quark potentials generated 
in the matching to PNRQCD enter this term as short- 
distance coefficients of the spatially non-local but tem¬ 
porally local, i.e. instantaneous, PNRQCD interactions. 
Ultrasoft interactions between the heavy quarks and the 
gluon field are contained in the last two terms of ([^ and 
the Aq terms in the first line. The latter do not contribute 
to colour-singlet production of a QQ pair, in which case 
ultrasoft effects enter for the first time at third order in 
non-relativistic perturbation theory. 

Perturbation theory in PNRQCD closely resembles 
quantum-mechanical perturbation theory, since the 
leading colour-Coulomb interaction is part of the un¬ 
perturbed theory. Thus, the propagator of PNRQCD in¬ 
cludes the leading Coulomb interaction exactly, which 
effects the required resummation of conventional per¬ 
turbation theory to all orders. The PNRQCD Feynman 
rules can again be found in jh). 

2.3. Multi-loop calculations 

The computation of matching coefficients with high 
precision involves fixed-order multi-loop calculations. 
For the topics reviewed in this article, the highest order 
is demanded in the calculation of the three-loop correc¬ 
tions to the static potential and the matching coefficient 
between QCD and NRQCD of the vector current. For 
calculations of this type it is necessary to have a high 
level of automation, which reaches from the generation 
of the Feynman diagrams to the final summation and ex¬ 
pansion in e. The scheme used for the three-loop vector- 
current matching coefficient requires particularly little 
manual interaction and is described in detail in ifT^ . 

For the generation of the contributing Feynman 
diagrams we use the Fortran program QGRAF IfTSlI . 
The output is passed via q2e lfT4l [TSll . which trans¬ 
forms Feynman diagrams into Feynman amplitudes, to 
exp lfT4l [TSll that generates FORM ifThl fTTl code. At 
the same time a topology is assigned to each diagram. 
The FORM code is processed to perform traces, apply 
projectors, and map the resulting scalar expressions to 
functions, which have a one-to-one relation to the to¬ 
pologies. They contain the exponents of the involved 
propagators as arguments. At this point one has in gen¬ 
eral a large number of different functions. Thus, in a 


next step one passes them to a program which imple¬ 
ments the Laporta algorithm lITSl and performs a re¬ 
duction to a small number of so-called master integrals. 
The latter have to be computed using analytical or nu¬ 
merical methods. In the case of the three-loop correc¬ 
tions to the static potential all but three coefficients in 
the e expansion could be computed analytically. For 
the three-loop corrections to the matching coefficient a 
significant number of integrals have been evaluated nu¬ 
merically using the program FIESTA mHH], which is 
based on sector decomposition. 

3. Top-quark pair production near threshold in e'^e~ 
annihilation 

Top-quark pair production near threshold in an¬ 
nihilation provides a unique possibility to measure a 
number of parameters with high precision. Among them 
are the top-quark mass and width, the strong coup¬ 
ling constant, and the top-quark Yukawa coupling. In 
particular the mass determination attracts lot of atten¬ 
tion, since the current uncertainty can be improved by 
about an order of magnitude, and at the same time the 
renormalization scheme is precisely defined. To reach 
this goal precise predictions are needed including both 
QCD and electroweak higher order corrections. The 
top-antitop system near threshold involves only scales 
in the perturbative regime Il22ll . which allows for sys¬ 
tematic approximations. Next-to-next-to-leading order 
(NNLO) QCD corrections have been computed at the 
end of the 1990s by several independent groups lfmi23l - 
l28l . The results are summarized and compared in Il29l . 
The inclusion of the NNLO corrections led to a signi¬ 
ficant shift of the cross section, and the theoretical un¬ 
certainty estimated through scale dependence remained 
larger than naively expected. The prospect of an ac¬ 
curate mass measurement at a future high-energy 
collider therefore makes it mandatory to push the theor¬ 
etical calculation to the third order. 

In this section we sketch the computation of the cross 
section cr{e^e^ tt-\-X) close to threshold to third order 
in non-relativistic perturbation theory. Normalized to 
the cross section for the production of pairs in the 
high-energy limit, ctq = 47ra^/(3s), it is given by 

cr(e+e“ tt + X) 

R = ^ 

O-Q 

= 127re^ Im = i H- fejj , (5) 

where e, = 2/3 is the top-quark electric charge, and the 
vector polarization function near threshold has the form 
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X G(E) + ... . ( 6 ) 


For simplicity of notation, we restrict this expression to 
the terms from a virtual photon, neglecting the Z-boson 
contribution. The full expressions are given in IQ. c,, 
and dv are NRQCD matching coefficients with perturb¬ 
ative expansions in as, which are defined through the 
expansion of the vector current in terms of non- 
relativistic fields given by 

/'’) '■ = (A ^ f'cr' D\ + .... (7) 

omj 

The terms proportional to £ = yfs-lnit in (j^ arise from 
expanding the prefactor Ijq^ - 1/s and from the l/m^ 
suppressed current in (|^, whose matrix element can be 
reduced to the one of the leading current by equation- 
of-motion relations. Thus, the central quantity in the 
non-relativistic representation is the two-point function 


G(£) = 




d‘^xe‘^-^° 


2NM- 

<0| T{ [«AV';r](0))|0)|NRQCD 


( 8 ) 


of non-relativistic currents, which has to be computed 
within NRQCD. 

In the following subsections we discuss various in¬ 
gredients which are needed for the third-order predic¬ 
tion to cr(e^e“ tt + X). In particular, we describe 
in Section 3.1 the calculation of the three-loop correc¬ 


tions to the matching coefficient Cv and in Section 3.2 


the three-loop corrections to the static potential. The 
latter is an important ingredient for the non-relativistic 
correlator which is discussed in Section 13.31 Ultrasoft 
corrections appear for the first time at third order and are 
discussed in Section 3.4 before presenting results for S- 
wave QCD contribution to the third-order cross section 
in Section |3.5| Section |3.6| contains further results, in 
particular a discussion of Higgs boson contributions at 
order aUs, of P-wave effects and of non-resonant con¬ 
tributions. 


3.1. Matching of the vector current, Cy 

One ingredient of the N^LO corrections to top-quark 
threshold production is the three-loop matching coeffi¬ 
cient between QCD and NRQCD. It is defined via the 
equation ESI 

Za.QCD Tqcd = Cv Z2,NRQCD Zj' TnrQCD , (9) 

where Z2 _qcd and Z2 ,nrqcd are the on-shell wave func¬ 
tion renormalization constants in QCD and NRQCD, re¬ 
spectively. Tqcd and Tnrqcd represent the amputated, 
bare electromagnetic current vertex functions evaluated 



Figure 1: Feynman diagrams contributing to C 3 . Straight and curly 
lines denote heavy quarks with mass m and gluons, respectively. 


for on-shell heavy quarks directly at threshold, that is, 
with zero relative momentum. It is understood that they 
are expressed in terms of the renormalized QCD coup¬ 
ling and the pole mass. One furthermore chooses for the 
photon momentum - 4m^. In this kinematic config¬ 
uration Tqcd corresponds to the hard part of the vertex 
function in full QCD. Moreover, we have Z 2 ,nrqcd = 1 
and only tree-level contributions to Fnrqcdj since loop 
corrections reduce to scaleless integrals in dimensional 
regularization. The renormalization constant Zj takes 
care of the infrared divergences of the renormalized on- 
shell QCD vertex leading to a finite result for c,,. 

One- and two-loop QCD corrections to c,, have been 
known since more than fifteen years and have been com¬ 
puted in on and ll^ 1^ [33l . respectively. The so- 
called renormalon contribution, consisting of the one- 
loop diagram with arbitrarily many massless quark loop 
insertions in the gluon propagator, has been computed 
in ES- More recently, the fermionic three-loop term 
became available Ha EH and the gluonic corrections 
have been computed in ESI- One-loop electroweak cor¬ 
rections in the SM are available from E3 and correc¬ 
tions of OiaUs) have been considered in Il38ll^ . Super- 
symmetric one-loop corrections to c,, have been com¬ 
puted in l40l . 

Sample Feynman diagrams contributing to Cy at the 
three-loop order are shown in Figure After generat¬ 
ing all relevant Feynman amplitudes and mapping them 
to scalar functions we use the program crusher mi 
for the reduction to about 100 master integrals. Some 
of them are known analytically ESiia but the ma¬ 
jor part is evaluated numerically using the program 
FIESTA lfT9H2TII . which provides for each coefficient 
in the e expansion an uncertainty originating from the 
Monte-Carlo integration. When summing the individual 
contributions these uncertainties are added in quadrat¬ 
ure. Such an approach requires strong checks which are 
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described in detail in ifT^I^ . A very powerful check is 
provided by the change of the master integral basis. We 
employ the integral tables generated during the reduc¬ 
tion procedure in order to re-express the master integ¬ 
rals, which are not known analytically, through differ¬ 
ent, in general more complicated ones. This transform¬ 
ation is done analytically for general space-time dimen¬ 
sion d. In a next step the new master integrals are again 
evaluated with FIESTA and inserted into the new expres¬ 
sion for c,,. For the final result we find excellent agree¬ 
ment within the uncertainties which is a strong check on 
the correctness of the result. 

In the following we present the result specifying Nc - 
3 and jj - m 


(«/) 

a\ 

c, ^ 1 -2.667— H- 
n 




[-44.551 -h0.407n,] 




^ J 



• 0.823 n 


- 2091(2)-H 120.66(0.0 l)n, 

H-singlet terms . ( 10 ) 


Results for individual colour coefficients and general 
values of p can be found in |[T2][35][36l. 


3.2. Coulomb potential, a 3 

In this subsection we summarize the potentials 
needed for top-quark pair production where our special 
emphasis lies on the three-loop corrections to the static 
potential. In momentum space the potential can be writ¬ 
ten in terms of operators expanded in 1/m. Its colour- 
singlet projection reads 


F(p,p') = -'Vcia,) 


47TCfas 


+ +0{l/m^), ( 11 ) 


m|q| 


where for brevity only the first two terms are shown; the 
1 /m^ operators can be found in lb). The coefficients of 
the operators are expanded in otj where to N^LO the po¬ 
tential coefficient ’Vdas) is needed to three-, 'V\/,„(as) 
to two-loop and the coefficient of the 1 /m^ term to one- 
loop accuracy. Note that there is no tree-level contribu¬ 
tion of order 1 /m^. 

The 1/m^ coefficients are known since long (see 0 
for a detailed discussion). Also ^i/m has been com¬ 
puted more then ten years ago Il43ll . however, only up to 
the constant term in the e expansion. Since the poten¬ 
tial loop-momentum integrals which have to be solved 
to obtain the corrections to the Green function G{E) 
are divergent, the 0{e) term of the two-loop correction 






Figure 2: Sample Feynman diagrams contributing to the static poten¬ 
tial at tree-level, one-, two- and three-loop order. Solid and curly lines 
represent quarks and gluons, respectively. In the case of closed loops 
the quarks are massless; the external quarks are heavy and treated in 
the static limit. 



Figure 3: Families of Feynman integrals needed for the calculation 
of « 3 . The solid lines correspond to relativistic massless propagators 
and the zigzag lines represents static propagators. 


to T^i/m is also required. It can be found in the Ap¬ 
pendix of Il44ll . The three-loop correction to 'Vq, com¬ 
monly denoted 03 , has been computed by two independ¬ 
ent groups msiii- The corresponding results shall be 
discussed in more detail in the following. 

In Figure sample Feynman diagrams contributing 
to 'Vc are shown. The external quarks are considered as 
heavy and thus for them the static approximation is ap¬ 
plied. The only momentum scale in the problem is given 
by the momentum transfer q between the quark and an¬ 
tiquark which has only non-vanishing space-like com¬ 
ponents. Furthermore, the static heavy-quark propag¬ 
ators have the form l/(-v ■ p - SjiO) with v ■ p - pQ 
and Sj = ±1 and thus do not depend on q. As a con¬ 
sequence, any Feynman integral that contributes to 03 
can be mapped to one of the three graphs shown in Fig- 
ure|^where solid lines stand for usual massless propag¬ 
ators of the form l/(-p^ - 10)“' and zigzag lines stand 
for static propagators which might also be raised to an 
integer index. In case the latter type of propagators 
is absent the integrals reduce to usual massless two- 
point functions, which can be treated with the help of 
MINCER l48l . Note, however, that the presence of the 
static lines significantly increases the complexity of the 
problem. 

In general the integrals involve up to fifteen propag¬ 
ators (including an irreducible numerator). In order to 
simplify the reduction problem we apply in a first step 
partial fraction identities to arrive at various subfamilies 
of integrals with at most three linear propagators. Thus, 
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Figure 4: Families with twelve indices (eleven propagators plus one 
scalar product) after applying partial fractioning to the “Mercedes” 
graph of Figure]^ 


any resulting integral is labelled by twelve indices one 
of which corresponds to an irreducible numerator and 
three indices correspond to linear propagators. As an 
example we show in Figure]^ the resulting families for 
the “Mercedes” graph of Figure]^ 

The next step is a reduction of all the Feynman integ¬ 
rals to master integrals by solving integration-by-parts 
relations na. This is done with the help of FIRE (Feyn¬ 
man Integral REduction) ||50H5^ . 

In total there are 41 master integrals which contrib¬ 
ute to the three-loop static potential. Ten integrals have 
no static lines and are thus known since long ||5TH55I| . 
Fourteen master integrals contain a massless one-loop 
insertion which can easily be integrated in terms of F 
functions using standard formulae. As a result one ob¬ 
tains two-loop integrals where one of the indices has 
a non-integer exponent involving the space-time para¬ 
meter e - {4- - d) 12. Results for these integrals are 
shown in |j56l. The 17 remaining integrals are more 
complicated and can only be computed in an expansion 
in e. Explicit results can be found in Ell (see also ED). 
Actually, all but three coefficients in the e expansion 
could be computed analytically. The three numerical 
results are known to a sufficiently high precision for all 
foreseeable phenomenological applications. 

An important check for the correctness of the calcu¬ 
lation is the use of a general QCD gauge parameter Hi. 
Note that the computational price one has to pay is quite 
high: a rough estimate of the complexity based on the 
number of integrals which have to be reduced to masters 
shows that the linear ^ term is about seven times and the 


term even 18 times more complicated than the Feyn¬ 
man gauge expression. Let us mention that nevertheless 
all occurring integrals could be reduced with the help of 
FIRE. 

We write the perturbative expansion of 'Vc in the 
form 


^c(as) 


H- 


/ ^ ' 

47r 

v3 

+ (^ 


\42r. 

1 ^ 


■d) 


4;r 




( 2 ) 


+ 0(at). 


( 12 ) 


Up to two-loop order the coefficients are finite, 
however, at three loops an infrared divergence oc¬ 
curs EH, which is related to ultraviolet divergences in 
the calculation of the ultrasoft corrections. It is con¬ 
venient to subtract these divergences in y(p, p') and add 
them back to the ultrasoft calculation. In the case of the 
Coulomb potential the infrared divergence is cancelled 
after adding the MS counterterm 




,C.t. 


asCp 3 as 
6e ’ 


(13) 


to y(p, p') in (11 1 . One finally arrives at 




( 0 ) _ 


c 

( 1 ) 


= 1 , 






( 2 ) _ 
C 

(3) 


2 2 

02 + (2ai/^o -i-^i)ln ^ -i-^gln^ 

= ^3 H- (2aif5i + P 2 + 3(32^0 + 


X In ^ -H H- 3 ai/?ojln^ ^ 

+ (14) 

q2 

Note the In term in W'® not associated with the QCD 
beta function is the remnant of the infrared divergence. 
We refrain from providing explicit results for a\ and 02 
which can be found in the literature ll43l l60H^ . even 
including higher order terms in e EHH. The result for 
02 reads isjiii 
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-56.83(1) 


^abcd ^abcd 
F F 

2NaTf 


«/■ 


502.24(1) C] + (-136.39(12)) 


^abcd ^abcd 

2NaTf 


dfcddabcd ^ ^4 _ ^^2 + 

INaTf ~ 96iv| 


^abcd ^abcd 


Nc{n^+6) 


2Na Tf 48 


( 15 ) 


( 16 ) 


As mentioned above, some of the ingredients entering 
03 are only known numerically; they enter four of the 
colour factors in ( [TSl l and are taken from ll46l l57l . Also 
parts of the N^LO Coulomb potential are known lIMI . 
but not needed for the third-order cross section calcula¬ 
tion. The two- and three-loop corrections to the static 
potential of heavy quarks in the colour-octet state have 
been computed in Il65l and ll66l . respectively. 

3.3. Third-order potential corrections to the PNRQCD 
correlation function 

Once the matching coefficients of NRQCD and PN¬ 
RQCD are determined, the correlation function G{E) 
defined in (|^ has to be calculated to the third order in 


PNRQCD perturbation theory. Since the leading-order 
Coulomb potential yo(Q in Q is part of the unperturbed 
Lagrangian, the position-space propagator is the matrix 
element Go(r, r';^) = {r|Go(£)|r') of the Green oper¬ 
ator Go(£) = [Ho-E-ie]^^. The propagator propagates 
a heavy quark-antiquark pair, and r refers to the separa¬ 
tion of the heavy quark-antiquark fields. Since one usu¬ 
ally works in the centre-of-mass (cms) frame, the cms 
motion is irrelevant. The unperturbed theory is not free, 
but it is still exactly solvable, since Hq is the Hamilto¬ 
nian of the Coulomb problem. An explicit expression 
for the Coulomb Green function in terms of Laguerre 
polynomials and a partial-wave decomposition is 

iiiiMl 


G 


( 1 , 8 ) 

0 


(r,r',E) 


'”^„-y(r+r') 

27r 


oo 


z 


(2/+l)(2yr)'(2yr')'P/(^j J 


s\Lf‘^^\2yr)Lf‘^^\2yr') 
(s21 + iy.(s-i-1 1-4)’ 


(17) 


where y - -\/-m{E -H ie) and A - (mas)l(2y) x {C/r (singlet); Cf - Ca/ 2 (octet)) for the propagation in a colour 
singlet and octet state, respectively. Further useful representations are summarized in IQ. The PNRQCD perturbation 


expansion of G(E) to third order reads 

G(E) = Go(£) -H 6iG{E) + SiGiE) h- d^GiE) h- ... (18) 

with Go(£) = <0|Go(F)|0) = G[|‘^(0, 0; E), and 

6,G(E) = <0|Go(£)ffiyi/Go(£)|0), (19) 

62 G(E) = <0|Go(£)ffiyi/Go(£)ffiyiiGo(£)|0)-H<0|Go(£)(5y2/Go(£)|0), (20) 

6iG{E) = <0|Go(F)ffiyi/Go(F)ffiyiiGo(F)ffiyi/Go(F)|0) -H2<0|Go(F)i^yiiGo(F)idy2iGo(F)|0) 

H- <0|Go(F)ffiy3iGo(F)|0) (5“G(£). (21) 


The corrections consist of potential contributions from 
single and multiple insertions of perturbation potentials 
6Vk of order G(y*) relative to the leading-order Cou¬ 
lomb potential, which arise from radiative corrections 
to the Coulomb potential and further non-Coulomb po¬ 


tentials such as ^i/m in (111. Starting from the third 


order, there is an ultrasoft contribution, 6‘‘^G{E), which 


is briefly discussed in Section 3.4 


The first term on the right-hand side of (21 1 repres¬ 


ents the triple insertion of the one-loop correction 6Vi to 
the Coulomb potential. The presence of four propagat¬ 
ors of the form turns this into a complicated expres- 
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sion. However, all the integrations are ultraviolet-finite, 
and this allows the triple insertion to be converted into 
sums that can be evaluated numerically as done in ||69|. 
Instead of computing order by order multiple insertions 
of 6Vi, one can also include 6V\ into the unperturbed 
Lagrangian and solve for the Green function of the as¬ 
sociated Schrodinger problem numerically, which sums 
the insertions of 6Vi to all orders. The comparison of 
the two approaches performed in ll69l shows that good 
convergence of the PNRQCD expansion to the exact 
result is only achieved when the renormalization scale 
ji is larger than 30 GeV. Somewhat surprisingly PNR¬ 
QCD perturbation theory already breaks down at scales 
larger than the natural scale of the inverse Bohr radius, 
but works when p is chosen closer to the hard scale. 
The explicit verification of this scale choice by com¬ 
parison with the exact result, which is only possible for 
the Coulomb potential, lends support to the scale choice 
that will be made in the evaluation of the full third-order 
cross section below. 

The second and third term on the right-hand side of 
(21 1 are associated with the single insertion of third- 
order potentials, such as the three-loop correction to the 
Coulomb potential discussed in the previous subsection. 


and the double insertion of a second-order potential to¬ 
gether with the one-loop correction to the Coulomb po¬ 
tential. While algebraically simpler than the triple inser¬ 
tion, they turn out to be more difficult to evaluate. The 
reason is that potentials more singular than 1 /r as r —> 0 
cause ultraviolet divergences of the potential integrals, 
which have to be consistently calculated and factorized 
in dimensional regularization, so that the divergences 
cancel with the matching coefficients and the finite term 
is correctly computed. At the same time, a calculation 
of the all-order diagrams summed by PNRQCD perturb¬ 
ation theory in d dimensions is not possible, since the 
propagator ( [T7] i is known only in four dimensions. 

The calculation of these non-Coulomb potential 
terms is described in detail in iTTOl . (The corresponding 
result for the bound-state parameters has already been 
presented in lITTI and for the Green function in ll72l . 
though no details were given in this work.) Here we dis¬ 
cuss as an example the single insertion of the 1/r^ po¬ 
tential to illustrate the strategy of the calculation. In mo¬ 
mentum space the tf-dimensional 1/r^ potential W'l/m is 
of the form where a is an integer. A 

single insertion is an integral 


I[x + ae] 



Go{pi,p2', E) 



—r] Go(P3>P4;£’), 


( 22 ) 


where q,y = p, - pj. Some of the potential insertions are multiplied by a divergent coefficient function w{e). This 
means that one should calculate I{x + ae\ to order e. However, these divergent coefficient functions always appear in 
conjunction with a counterterm with a slightly different momentum dependence, such that the potential expanded in e 
is finite. We therefore consider the expression 


1 


(q2)^ 


w(e)- 


with 


w{e) = 


w 


,(l/e) 


+ w + + 0{e^), 


(23) 


and define the corresponding counterterm-including single-insertion function as 


1 


J[x + ae-,w(e)] - -w^^^^^\I[x + ae] - I[x]\ + \w+ w'^^^e\I{x + ae\. 




(24) 


The advantage of this expression is that it avoids the 
need to calculate in I{x + ae] the a-independent 0{e) 
terms, since they drop out in the difference in brackets in 
the hrst term. These terms would indeed be diificult to 
obtain. The second term is multiplied by a hnite series, 
hence the 0{e) term of I{x] is never required to obtain 
the hnite part of 7 as e ^ 0. 

The single insertion of the 1/r^ potential generates 


ultraviolet 1 je poles from the integration over the po¬ 
tential loop momenta, which are related to the singular¬ 
ities in the dimensionally regulated hard current match¬ 
ing coefhcients, and which have to be properly factor¬ 
ized. Power counting shows that the insertion has an 
overall divergence coming from diagrams with less then 
two gluon exchanges, and a vertex subdivergence, when 
there is no gluon exchange between the external vertex 













M. Beneke and M. Steinhauser, Non-relativistic high-energy physics: top production and dark matter annihilation (2015) l-\36\ 


10 


and the potential insertion. To accomplish the correct parts (and show the corresponding diagrams below): 
factorization, we divide the integral into four different 


/[- + ae] - /fl[- + ae] + 



Similar notation applies to the counterterm-including 
insertion function J. The first two terms map to two- 
and three-loop diagrams, which carry an overall di¬ 
vergence and a divergence in the vertex subgraph(s) 
without gluon exchanges. The fourth term represents 
a sum of diagrams to all orders, but is finite. The third 
term has a divergence in the left vertex subgraph and all¬ 
order summation to the right of the potential insertion. 


where the notation “ > 1” refers to all ladder diagrams 
summed by Go(E) with more than one gluon exchange. 

The calculation of the first two parts is straightfor¬ 
ward, since it involves only ordinary, solvable, dimen¬ 
sionally regularized two- and three-loop diagrams. Nev¬ 
ertheless, the divergent part must be properly factor¬ 
ized. To be explicit, the two-loop integral for diagram 
(a) evaluates to 


+«] = (r)' 


y-mE I mE\"l mE\ 
\ / I 4.7r / 


G[,°“^(Pi,P2;£)G'°“>(P3,P4;£) 


,2 

-2" r(l - e - M)r(e M)2r(-i +u + 2 e) 


(qly 


r(i-e)T(2e + 2u) 


(26) 


from which 7a[l/2 + ae;w(e)] follows according to the definition (24 1 . The overall y-niE factor indicates that the 
divergence of this integral persists in the imaginary part of the correlation function and is proportional to Gq*'“^(£’). 
Since the divergent part of the hard matching coefficient Cy multiplies the (f-dimensional correlation function, we must 
write the pole part of Ja[\ +ae', w(e)] in such a way that it multiplies the li-dimensional expression for G^q‘^'’‘\E), which 
is given by (jf - j {An)) 




■f 


d‘‘ *p 


-1 


[ yj-mE 


{2n) 


d-l 


87r3/2 


mE 

1,2 


e^r.r(-l + e). 


(27) 


This results in (L^ = - ^ ln(-4m£’/p2)) 




+ae-, w{e)] 


2n^(a + l)e2 ° 
m^Cfas 




•j(2fl(ln2-l)w('/^^-Hwj 


2n^{a + l)e 




8n^(a + 1)4 


(i/f) 


ln^2-21n2- 


24 


(2a + 5) -2a - 2{a - 1 - l)Lx 


-{a+ f)E\ I -I- wl ln2 - (fl -I- 1)L2 - 2 - a I-w' 




(28) 
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where the 1/e pole parts now multiply the d- 
dimensional zero-exchange Coulomb Green function. 

Contributions such as part c, which have subdiver¬ 
gences and represent all-order graphs are the most com¬ 
plicated ones, since that appears to the right 

of the insertion is not known in (f-dimensions. The 
strategy consists of isolating the subdivergence at the 
integrand level, so that it can be factorized without ever 
requiring an explicit representation of The 

hnite remainder can then be evaluated in four dimen¬ 
sions. This strategy always works, because the diver¬ 
gences can only arise from subdiagrams with a hnite 
number of loops. This must be so, since the divergences 


must cancel with infrared divergences in the matching 
coefficients, which are computed in hxed-order perturb¬ 
ation theory. On the technical level, the hniteness of 
ladder diagrams with potential insertions once there are 
sufficiently many ladder rungs follows from the fact that 
every rung reduces the degree of divergence. For ex¬ 
plicit results for the other parts of 7[ | -H ae; w(e)] we 
refer to uni. An important consistency check is that the 
factored divergent parts of Ja, 2Jt and 2Jc are all the 
same. This is necessary for the sum of all contributions 
to add to a term proportional to the full Green function 
Go(£) = Gf^\E) + g[)‘“^(£) -h g[,^‘“’(£). That is, for 
the sum of all parts we have 


, (i/e) m(2a(ln2 - l)w^'^''^-H w) 

J[-+ae-w{e)] = n 2 ^ 0 ^-^) + —-- ’-Go(E) + 0{e\ (29) 

2 2n-^(a + l)e^ 27T‘^{a -H l)e 


as is required for cancelling the divergent part with 
the hard matching coefficient multiplying Gq{E). The 
method also applies to double insertions, where some 
expressions require the factorization of a 1 /e pole mul¬ 
tiplying the single insertion of the NLO Coulomb poten¬ 
tial. In the end, one verihes the cancellation of all poles 
in the sum of all terms, and evaluates the remainder nu¬ 
merically. 

3.4. Ultrasoft correction 

A contribution from the ultrasoft loop momentum re¬ 
gion has to be taken into account for the hrst time at 
the third order in perturbation theory. The ultrasoft cor¬ 
rection 6‘‘^G{E) to pT] ) is technically and conceptually 
the most complicated contribution to the third-order cor¬ 
rection to G(E), since its (sub) divergence structure is 
rather involved and the integrals are more difficult. The 
correction to the Green function relevant to the top-pair 
cross section was computed in ITTTI . with the subtrac¬ 
tion scheme already developed in ca to determine the 
bound-state residues. In this subsection we want to sum¬ 
marize the most important features. 


The ultrasoft interaction terms in the PNRQCD Lag- 
rangian 0 relevant to third order are given by 

gsil/\x)\AQ{t, 0) - X ■ E(f, 0)]i/f(v) 

+ g,x\x)[Aa{U 0) - X ■ E(f, 0)];r(^)- (30) 

The leading Ao(f, 0) couplings can be removed by a 
redehnition of the heavy-quark helds with a time-like 
Wilson line, which modihes the production current. The 
Wilson lines cancel for colour-singlet currents, hence 
the Ao(f, 0) couplings have no effect on top-pair produc¬ 
tion in collisions. Note that x ~ 1 /v and g^E ~ 
for ultrasoft gluon helds and thus the remaining chro¬ 
moelectric dipole interaction is suppressed by rel¬ 
ative to the kinetic term in the action. Two ultrasoft in¬ 
teraction vertices are required to form an ultrasoft loop, 
from which it follows that the leading ultrasoft contri¬ 
bution arises only at the third order. 

Employing the Feynman rule for the x ■ E(f, 0) ver¬ 
tex, the ultrasoft correction to the Green function can be 
expressed in the form 


d“G(£) = ig]C, 




d^k 

( 2 ^ 


kg r ■ r' - (r ■ k)(r' ■ k) 
+ ie 


G;/>(0, r; E)Gf \r, r'; E - ko)G«(r', 0; E) 


, (31) 
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with the understanding that one picks up only the pole 
at k** = |k| - ie in the gluon propagator. However, 
this expression cannot be used in practice, because it 
is ultraviolet divergent. Instead, one reverts the deriv¬ 
ation of the chromoelectric dipole interaction and uses 
the Feynman rules for the threshold-expanded NRQCD 
momentum-space diagrams together with dimensional 
regularization. 

The ultrasoft correction 5‘‘^G(E) has ultraviolet diver¬ 
gences from the integral over the three-momentum k of 


the ultrasoft gluon, and from the subsequent potential 
loop integrations. The former take the form of a single 
insertion of a third-order potential and of a one-loop cor¬ 
rection to the coefficient dy of the 0{v^) suppressed de¬ 
rivative current. Note that ultraviolet divergence has to 
be regularized in dimensional regularization to be con¬ 
sistent with the calculation of potential insertions and 
hard-matching coefficients. For this reason it is conveni¬ 
ent to introduce appropriate subtraction terms to cancel 
the ultrasoft subdivergences, which leads to EEl 




rp2.]2 

r d‘^-^e d‘’-^e' f 

[R \ 

J (27r)"'-i (271^-^ 1 


(- 1 )- 


f + t'^ 


6m?- 




['^ 17 : 


{2nY-^ {2nY-'^ 




6U - 6Vc,. 


iG^l\^'J'-E) 


(32) 


where SVc.t. represents the potential subtraction, and 
SU is the ultrasoft insertion (containing the octet Green 
function). The first line of ( |32] i is related to the renor¬ 
malization of the G(v^) suppressed vector current. In 
fact, ddf'' contains the infrared divergence that was sub¬ 
tracted to obtain the finite expression for dv in (|^. The 
divergences in the second line are associated with po¬ 
tential insertions. The counterterm ( [T3| ) that is needed 
to make the three-loop Coulomb potential finite after 
coupling renormalization is contained in dVc.i. above. 
The remaining divergences of (32i from the potential 
loop integrations are associated with the three-loop hard 
matching coefficient Cy and isolated as explained in ll74ll . 
The final result is then evaluated by a combination of 
analytical and numerical methods. 


3.5. Third-order cross section 

In this subsection we put all third-order corrections 
together and briefly discuss their numerical effect on 
the production cross section of top-quark pairs (see 
also CS). We restrict ourselves to pure QCD correc¬ 
tions and furthermore only consider S-wave contribu¬ 
tions. 

The (normalized) total cross section is given by 0. 
We parametrize R in terms of the potential subtracted 
(PS) mass lESl. This avoids the infrared sensitivity of 
the pole mass CiIItHI, which would otherwise prohibit 


a mass determination with accuracy below Aqcd- The 
PS mass is related to the pole mass via 

= m, - 6m,(jif), (33) 

where 6m,{fif) is obtained from the static potential solv¬ 
ing the following integral 

6m,(i,f)^-\ f 7^V(q). (34) 

The factorization scale /Jf is part of the definition of the 
PS mass. In the following we set it to /r/ = 20 GeV. To 
obtain 6m,{fif) to N^LO ll6^ the three-loop coefficient 
of the static potential (cf. Section [3^ is needed. 

We are now in the position to numerically evaluate 0 
where our default input parameters are given by 

aAMz) = 0.1184, 

= 171.3 GeV, 

F, = 1.4 GeV. (35) 

As the central value for the renormalization scale we 
adopt p = 80 GeV. 

In Figure the normalized total cross section is 
shown as a function of -sfs in the threshold region. The 
bands are obtained by simultaneous variation of the 
renormalization and factorization scale between 50 and 
350 GeV. After the inclusion of the third-order correc¬ 
tions one observes a dramatic stabilization of the per¬ 
turbative prediction, in particular in and below the peak 
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VT (GeV) 

Figure 5: R as a function of the e'^e~ centre-of-mass energy in the 
threshold region. LO, NLO, NNLO and N^LO results are shown in 
grey, blue (lowest band at ^/s = 347 GeV), red and black where the 
bands indicate the uncertainty due to scale variation. 

1.25 
1.20 
1.15 
^ 1.10 
1.05 
1.00 
0.95 
0.90 

343.6 343.8 344.0 344.2 

(GeV) 

Figure 6: LO, NLO, NNLO and N^LO results for the height (vertical 
axis) and position (horizontal axis) of the peak. The error bands ori¬ 
ginate from scale variation and uncertainties in a,. 

region, where most of the sensitivity to the top-quark 
mass comes from. In fact, the N^LO band is entirely 
contained within the NNLO one in this region. This is 
different above the peak position where a sizable neg¬ 
ative correction is observed when going from NNLO to 
N^LO. For example, 2 GeV above the peak this amounts 
to about -10%. 

It is interesting to have a closer look at the position 
and height of the peak. Figure shows the peak height 
and its position at LO, NLO, NNLO and N^LO where 
the error bars reflects the uncertainty due to the scale 
and as variation, which are added in quadrature. In 
the position of the peak one observes a relatively big 



(GeV) 

Figure 7: N^LO corrections to R(s) normalized to its value for = 
80 GeV where the band corresponds to scale variation between 50 
and 350 GeV. For the (red) solid curve the input value for m)"® in j35| 
by -h 50 MeV. 


jump from LO to NLO of about 400 MeV, and from 
NLO to NNLO of approximately 150 MeV which re¬ 
duces to only 50 MeV from NNLO to N^LO. Further¬ 
more, the NNLO and N^LO uncertainty bands show a 
significant overlap. The uncertainty of the peak posi¬ 
tion amounts to about 70 MeV, 60 MeV and 30 MeV 
at NLO, NNLO and N^LO. As far as the height is con¬ 
cerned there are big jumps from LO to NLO and NLO 
to NNLO, however, the N^LO correction stabilizes the 
perturbative series and a shift below 3% is observed. 

In Figure]^ we study the sensitivity of the total cross 
section on the variation of the top-quark mass. This is 
done by normalizing the N^LO band of Figure |^to the 
result obtained for /r = 80 GeV. The hatched band in 
Figure represents the theoretical uncertainty of this 
normalized quantity. The red curve in Figure is ob¬ 
tained by adopting = 80 GeV and shifting by 
-1-50 MeV. Above threshold basically no effect is ob¬ 
served. However, the strong variation of the cross sec¬ 
tion below threshold implies a several percent effect 
even for such a small variation of m)’®, and provides a 
handle for the precise determination of the top-quark 
mass. From Figure |7] one may expect a precision of 
about 50 MeV, and it is obvious that the small theor¬ 
etical uncertainty after including the third-order correc¬ 
tions to R is crucial to reach this goal. Similar results 
are obtained for variations of other parameters like the 
strong coupling constant and the top-quark width ll75l . 

Finally we mention that (almost complete) next- 
to-next-to-leading logarithmic predictions for the total 
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cross section have been obtained GSlIlol. In this ap¬ 
proach, the logarithmically enhanced terms in v of the 
third-order correction are included, but not the sizable 
“constant” terms. On the other hand, logarithms of v 
are summed to all orders. At first glance there appears 
to be good agreement with our N^LO results. However, 
a detailed comparison of the different terms included in 
the two approaches is required to determine whether the 
agreement is more than coincidental. 

3.6. Further results 

With the third-order QCD corrections to the dominant 
S-wave production mode completed and uncertainties 
under good control, attention must be paid to other cor¬ 
rections that may be less challenging from the technical 
perspective, but important phenomenologically. There 
are QED and electroweak corrections, where we count 
tTem ~ ctew ~ ~ v^', Higgs boson contributions; pro¬ 

duction of tt in a P-wave state; non-resonant contribu¬ 
tions to the physical final state W^W^bb', and a con¬ 
sistent treatment of initial-state radiation with controlled 
dependence on the factorization scheme for the initial- 
state electron distribution function. We discuss some of 
these items below (see MM for a general discussion 
of initial-state radiation in this context and some relev¬ 
ant formalism). 

3.6.1. P-wave contribution 

The P-wave contribution arises from the axial-vector 
coupling of the Z-boson to the tt pair. Since the ima¬ 
ginary part of the axial-vector polarization function 
n(“)(gr2) is 0{v^) suppressed in the non-relativistic limit, 
a NLO calculation of in PNRQCD results in 

a third-order contribution to the total top-pair produc¬ 
tion cross section near threshold. The NLO calculation 
was performed in 18^ The NLO correction stabilizes 
the scale dependence, but the result is ambiguous until 
a scheme is defined to combine the resonant with non¬ 
resonant contributions, see Section [3. 6. 3| below. How¬ 
ever, the P-wave contribution is always a very small 
contribution, below 1%, to the total cross section. 

3.6.2. Higgs contributions 

The Higgs contributions are of particular interest, 
since they provide sensitivity of the cross section 
cr(e+e“ ^ it + X) to the top-quark Yukawa coupling 


'Some results for the P-wave Green function were already ob¬ 
tained in I28ll84tl8^ . but none of these computations were performed 
in dimensional regularization as required for consistency with other 
pieces of the calculation. 


to the Higgs boson. The dependence originates from 
Higgs exchange between the produced top quarks (po¬ 
tential region), and from radiative corrections to the pro¬ 
duction vertex (hard region). These effects are incor¬ 
porated into the PNRQCD framework via modifications 
of the vector current which result from simultaneously 
integrating out the top quark and the Higgs boson, as¬ 
suming Mh ~ mt for power counting. Lurthermore, 
a new operator occurs in the effective theory, which in 
momentum space is given by HSl 




amnf 

y 2 ^2 11^2 ’ 


(36) 


where siy is the sine of the weak mixing angle and 
is the W boson mass. In coordinate space, this is equi¬ 
valent to a delta function potential that approximates the 
short-range Yukawa potential. If we employ the count¬ 
ing rule a ~ it is easy to see that 6X.h gives contribu¬ 
tions which are parametrically of the same order as the 
ones from the third-order PNRQCD Lagrangian. 

We parametrize electroweak corrections to the match¬ 
ing coefficient Cy as 




Cv = 1 + —rC + 


77-2 c2 


(37) 


where the the ellipses stand for QCD corrections which 
can be found in ( [T0| ). The complete one-loop elec¬ 
troweak corrections to the matching coefficient c™ have 
been computed in llJTl . The Higgs boson contribution 
can be cast in the form 




3y 


(-2 + 5yl 


2-9yj,+ Uy- 2 
-Iny^ 


48y 


H 


•6y4) 




yi 


arctan 




(38) 


where yn - mtlMn and m, is the on-shell top-quark 
mass. Mixed contributions of order aas involving the 
Higgs boson have been computed in lf38ll where various 
expansions have been considered. Lor ^125 GeV 
an expansion around m, ^ leads to the best results 
for which can be written as 


//,mix _ 
^\\lh ~ 


rUf 


2 2 

(1 -yH,ii)ln — 

8 /i^ 


5.760 


+ 5.533yH,i/, - O.niy^i, + 0.0124y^i, 

+ 0.0304y^i, + 0.0296y^i, + ... , (39) 
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/i" (GeV) 

Figure 8: N^LO corrections to R{s) normalized to its value for p = 
80 GeV where the band corresponds to scale variation between 50 and 
350 GeV. For the dashed curves the input value for the top Yukawa 
coupling in the ttH vertex is rescaled relative to the SM value by a 
factor 1.5, 1.2, 0.5 (from top to bottom) while keeping the top-quark 
mass at its input value ID The central value to which the variation 
is normalized includes the SM Higgs contributions. 

withyH.ii = (1 -y]i)- 

The PNRQCD computation of the Higgs contribu¬ 
tions to cr(e+e“ —> tt + X) has been performed in |[87l 
and results in corrections at the 10% level as expected 
from the discussion of the corrections to the bound-state 
residue Il38l . Higgs contributions of order aa, clearly 
have to be incorporated in order to reach a precision at 
the few-percent level. 

Figurej^shows the effect of a modification of the top- 
Yukawa coupling relative to its SM value of the top-pair 
production cross section near threshold in a format sim¬ 
ilar to Figure]^ Comparing the width of the uncertainty 
band to the shift of the curves provides a first estimate of 
how well the top-Yukawa coupling can be constrained 
at a future e^e~ collider before the cms energy for open 
ttH production is reached ISTl . 

3.6.3. Non-resonant contributions 

The most important piece that is still missing for a 
realistic prediction of the total cross section is related to 
the fact that the top quark decays rapidly. The measure¬ 
ment of top-pair production near threshold really refers 
to the measurement of —> W^W^bb with cms en¬ 

ergy s/s near 4mj. The pure QCD contribution, which 
we referred to above, is definedhy the polarization func¬ 
tion evaluated in the complex energy plane with the pre¬ 
scription E ^ E ->r iYt. 

The limitations of this approximation manifest them¬ 
selves within the (NR)QCD calculation itself. The cur¬ 
rent correlation function GiE) exhibits an uncancelled 


ultraviolet divergence from an overall divergence of the 
form [5G(£’)]overaU OsEle in dimensional regulariz¬ 
ation Ga. Since E acquires an imaginary part F, ~ 
through the above prescription, the divergence 
survives in the cross section, 

(1'cQ'fw 

Im [dG(£')]overall X -. (40) 

e 

The divergence appears first at NNLO (since at LO, 
G{E) ~ V ~ Os), and results in dependence on an addi¬ 
tional regularization scale after subtraction. A con¬ 
sistent calculation therefore requires that one considers 
the process —» W^W^bb within unstable-particle 

effective theory Il88l[89l including the effects of off-shell 
top quarks and processes that produce the W^W^bb fi¬ 
nal state with no or only one intermediate top-quark 
line. The physical cross section is then the sum of two 
terms, 

CTe+e-^W+W-hb - 0'e*e-^[tt],cs(Ew) 
pure (NR)QCD 

Both terms separately have a “finite-width scale depend¬ 
ence”, and only the sum is well-defined. The calcula¬ 
tion of the non-resonant part has been performed so far 
only at NLO Il90l 1^ . including the possibility to ap¬ 
ply invariant-mass cuts on the top decay products. Parts 
relevant to NNLO are known ll9TH94ll . but the complete 
computation is still missing. For a further discussion of 
non-resonant contributions and unstable-particle effect¬ 
ive theory we refer to the review 1951 in this volume. 

4. Threshold resummation in hadronic collisions 

The production of top-quark pairs, and pairs of heavy 
particles in general, in proton-proton (LHC) or proton- 
antiproton (Tevatron) collisions is different from 
annihilation in two important aspects. First, the centre- 
of-mass energy of the colliding partons is not fixed. 
Non-relativistic theory is strictly relevant only in the 
kinematic region of small invariant mass of the heavy- 
particle pair near the kinematic limit. Second, the 
initial-state particles are coloured, so that the heavy- 
particle pair can be produced in different colour states. 
This implies that soft gluon effects are much more im¬ 
portant than in annihilation, both due to initial- 
state radiation, and due to the coupling to the coloured 
final state. (Following standard terminology, what we 
call “soft” in this section, refers to the ultrasoft region 
in previous sections.) In this section we discuss the the¬ 
oretical formalism that brings together non-relativistic 
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theory with the theory of soft-gluon resummation, the 
latter having been developed first for Drell-Yan pro¬ 
duction ll96l l97l and then extended to di-jets or pairs 
of heavy particles ||98U104I . We then consider the 
invariant-mass distribution of hadronically produced top 
pairs, and the total cross section for top and superpartner 
particle pairs. 

4.1. Joint soft and Coulomb resummation 

The theoretical quantities of interest are the hard- 
scattering cross sections & for the partonic subprocesses 

p(ki)p'(k2) H(pi)H'(p2) -I- X (42) 

with pp' e {qq,qq,qq, gg, gq, gq}, when the partonic 
centre-of-mass energy s is close to 4M^ = (mn -i- mn’)^, 
such that /3 - {1 - Isf-I^ is small. The two heavy 
particles H, H' can be in arbitrary, not necessarily 
identical colour representations R, R'. 

The expansion of & contains enhanced terms of the 
form ln^/3)", primarily (but not exclusively) related 
to soft gluon effects, and the “Coulomb singularities” 
(aj/yS)", which both cause a breakdown of the perturb¬ 
ation expansion for small [3. Although the power-like 
Coulomb singularities are formally stronger than the 
logarithmic soft-gluon effects, the importance of one or 


the other effect depends on the initial state (quark vs. 
gluon) and colour of the final state (for instance, sing¬ 
let vs. octet). Joint soft and Coulomb resummation 
0105111061 concerns the question whether both effects 
can be resummed simultaneously. 

4.1.1. Factorization 

This is a non-trivial issue, since the energy of soft 
gluons is of the same order as the kinetic energy Mff' 
of the heavy particles produced. Soft-gluon lines may 
therefore connect without parametric suppression to the 
heavy-particle propagators in between the Coulomb lad¬ 
der rungs, as well as to the Coulomb gluons itself (by 
virtue of the gluon self-coupling), impeding the stand¬ 
ard factorization arguments that assume either no coup¬ 
lings to the final state, or to the energetic initial-state 
particles. 

Factorization was investigated in 0105111061 . where 
it was shown that the partonic cross section factor¬ 
izes into three separate contributions, related to hard, 
soft and Coulomb effects. The latter two are coupled 
only through a convolution in an energy variable u, 
which roughly speaking accounts for the fact that near 
threshold the production of the heavy particle pair is 
sensitive to the small amount of energy radiated into soft 
gluons. The factorization formula reads 


i Ra 


Here E - s/S- 2m, is the energy relative to the produc¬ 
tion threshold. Referring to lfT06l for the derivation, we 
explain here the elements of the formula and properties 
of the effective interactions that lead to this simple form. 
The formula as written applies to the total cross section. 
A similar result without convolution integral holds for 
the distribution in the HH' invariant mass. 

The first factor accounts for the short- 

distance production of HH'. It depends on the masses 
of the heavy particles, the partonic initial state, but not 
on the small scales E and Mfi in the problem. The su¬ 
perscript i refers to a colour decomposition to be dis¬ 
cussed below. The function Jr„{E) accounts for the 
non-relativistic effects from the potential region and co¬ 
incides with the zero-distance Coulomb Green function 
G{E), generalized to the case of unequal-mass particles 
and an arbitrary irreducible SU(3) representation. At 
leading order, this is easily accomplished by replacing 


jRjE-‘^)Wf"(cj,p). 


(43) 


m by the reduced mass, and the colour-singlet coeffi¬ 
cient Cf of the Coulomb potential by is de¬ 

termined by 


^(R)brr(R')b pRa _ r\ pRn 

^a2C2 ^ciC 2 a 3 a 4 ''^^1^101010304)’ 


(44) 


where P^° projects on the irreducible representation R^ 
in the tensor product R ® R', and are the SU(3) 
generators in representation R. Finally, W^°{io,p) is the 
soft-gluon function in the representation 


Eq. (43 I holds up to corrections 0(J3^) for the total 
cross section near threshold. Within this approxima¬ 
tion the interactions of soft gluons with collinear and 
non-relativistic particles take a particularly simple form, 
which is essential for proving factorization. The coup¬ 
ling of soft gluons to an energetic initial-state quark with 
four-momentum in the light-like n direction is given by 
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the term 

^c(x) igsn ■ A,(;c+) | (45) 

in the position-space SCET Lagrangian, which is equi¬ 
valent to the eikonal approximation. Here denotes the 
n-collinear quark held, and h is another light-like vec¬ 
tor with n ■ h - 2. It is important that only the n ■ 
component appears, and the soft held is evaluated at the 
point x^l - {h ■ xl2)rf = x^rf due to the light-cone 
multipole expansion 0107111081 . This allows the soft- 
gluon coupling to be removed from the Lagrangian by 
performing the held redehnition ^dx) - Sn \x-)^f\x) 
IIIo3, where 


SfHx) = Pexp 



0 

dtn ■ A^(x + nt)T 


oo 


(R)c 


(46) 


denotes a light-like, soft Wilson line for a particle in 
the representation R of SU(3). Similar results apply to 
initial-state gluons. The leading soft-gluon interaction 
with the non-relativistic helds for H and H', 


id\x) igA°(xo) Hx) + t/'' (x) igsAs(xo) ij/'ix), (47) 


involves only the Aq helds at point (xq, 0), and therefore 
can be removed by i^/ix) = 5!if\xo)iA**’'(x), where w - 
(1,0) is a time-like four-vector, and 


5f'(x) = Pexp 


/-»oo 

igs I dtw A^x + wt) 
Jo 


(48) 




is the corresponding Wilson line. The decoupling of 
soft gluons from the Coulomb gluons in ladder dia¬ 
grams now follows from the identity = 

which in turn implies [i/r'''T*^^°i/fj(x -H r) = 

5®^^(xo) -I- r). Since the Wilson line 

in the adjoint representation is real and hermitian, and 
the argument is xq independent of r, the Wilson line 
drops out from the Coulomb interaction 

J iV [i/r'T(«)>](x + r) y [tA' ' T(*'^>'](x), (49) 

when expressed in terms of the redehned helds 01061 

[m. 

The soft gluons do not disappear completely. The 
short-distance production amplitude of the heavy- 
particle pair in the hard 2 —> 2 process is described by 
an operator of the form 

which is local (up to collinear Wilson lines Wc in collin- 
ear helds such as (^c for the initial-state partons). After 
squaring the amplitude and introducing the redehned 
helds, eight soft Wilson lines are left over. Since the 
effective Lagrangian, expressed in the redehned helds, 
does not contain soft interactions with the other helds at 
leading power, these Wilson lines can be collected into 
the soft function 




(51) 


r 


where T (T) denote (anti) time-ordering. Latin indices 
refer to SU(3) colour and the Wilson lines are under¬ 
stood to refer to the representations of the correspond¬ 
ing particles. Curly brackets (a) denote a colour multi¬ 
index, here aia 2 a^a 4 . 


is produced at rest at threshold, suggests that the gluon 
coupling to the pair can be described by a single soft 
function. Indeed, by virtue of 


o (tt) c{R') _ 

'-00102^ v.aiii V,fl2*2 ~ v,ap'~'pbibi’ 


:(R') 


(52) 


Lurther simplihcations that lead to (431 are closely 
connected to properties of the soft function 110511 . The 
physical intuition that soft gluons should couple only 
to the total colour charge when the heavy-particle pair 


where denotes the Clebsch-Gordan coefficient 

that couples the representations R, R' to the irreducible 


representation Ra inR^R', the soft function (51 1 can be 
related to a sum of simpler soft functions of the form 




R,A 

v,Ka 


]( 0 )| 0 ). 


(53) 
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Ms 


Figure 9: Sketch of the resummation of soft gluon corrections using 
renormalization group equations. Figure from (Ml. 


The Clebsch-Gordan coefficients can also be used to 
construct an orthonormal basis of colour tensors c®j for 
the hard-scattering amplitude with initial-state partons 
in the representation r ® r' = and HH' in R®R' - 
Za Ra, such that Jlpp'{a} = Zi ^hc index i enu¬ 

merates all pairs P, = {rp,Ra) of equivalent represent¬ 
ations rp and Ra that occur in the above decomposition 
of the tensor-product representations. For example, in 
top-pair production through gluon-gluon fusion, there 
are three combinations P; e {(1,1), (8s,8), (8 a,8)). 
Colour conservation implies that 



1 

yAimirp) 


p~'Ra* 


form an orthonormal basis. Defining 


(54) 


= M) pK 
[aa] ~ 


(55) 


the soft function (531 can be represented as a matrix in 
this basis: 


The formula now follows from the fact that 

can be shown to be a diagonal matrix to all 
orders in perturbation theory for all relevant cases M105L 
The factorization formula ( |43] l justifies earlier treat¬ 
ments of threshold resummation for heavy particles, 
where soft-Coulomb factorization has been put in as 


an assumption miiHmi, or the Coulomb-enhanced 
terms were not summed and technically considered 
as part of the hard function CMI. Soft-gluon re¬ 
summation in heavy-particle production has been per¬ 
formed in Mellin-moment space in previous works. We 
therefore note that the convolution in •ED turns into 
multiplicative soft-Coulomb factorization &pp'(N,iJ.) 

2 / H‘^p, JrA^) of the partonic cross sec¬ 

tions in moment space. 

Away from the production threshold large soft-gluon 
logarithms can appear in various other kinematic re¬ 
gions. When the relative velocity of the heavy particles 
is not small. Coulomb effects are not enhanced, and fac¬ 
torization takes the more standard form H ■ W without 
the J function. However, the hard and soft functions 
now depend on the kinematic invariants of the general 
2^2 scattering process, as do the anomalous di¬ 
mensions relevant to resummation. The soft function 
is no longer diagonal. Near threshold, interactions that 
change the colour-state of the heavy-particle pair appear 
only at the level of 00^) (0(J3)) corrections to the total 
cross section (amplitude). Some of these effects are dis¬ 
cussed in II106L but no calculation of 0(J3^) corrections 
in non-relativistic perturbation theory has so far been 
performed for hadronic processes. 

4.1.2. Resummation 

In the present framework the resummation of “Cou¬ 
lomb singularities” (aj/l)" is automatic in the J func¬ 
tion. The logarithms (as 10 jB)" are summed by solving 
the renormalization group equations to evolve the hard 
functions N from the scale jih ~ M, and the soft func¬ 
tions W from the scale ^s ~ M/0 to a common scale n /, 
chosen to be the factorization scale of the parton distri¬ 
butions, as illustrated in Figure]^ 

For the systematics of the combined resummations of 
the two types of corrections we count both Oslnfl and 
aj/jS as quantities of order one and introduce a paramet¬ 
ric representation of the expansion of the cross section 
of the form 


Cpp' 



lnj3go(as In/i) -i- gi(a* In^S) H- asgiCa, In/?) -H ... 
(LL) (NLL) (NNLL) 


X { 1 (LL,NLL); as,j3 (NNLL);...). 


(57) 
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The resummed cross section at LL accuracy includes 
all terms of order 1 /yS^ x a"^*ln^"y 6 relative to the 
Born cross section near threshold. Next-to-leading 
summation includes in addition all terms of order 
asln/3; a^{l//3 x In/?,In^jS);..while furthermore all 
terms 0 -^; a^{l//3, In^'* /?);... are included in the NNLL 
approximation. With this counting, the formalism de¬ 
scribed above is limited to NNLL, since 0(p^) terms are 
required beyond this order. A caveat applying to current 
NNLL treatments should be mentioned. Higher-order 
non-relativistic potentials cause ultraviolet singularities 
in the Coulomb function J, which cancel with the hard 


functions, and cause non-relativistic logarithms begin¬ 
ning with \n/3 at NNLL order. Present NNLL results 
(as described below) sum soft logarithms to all orders 
by evolution of the soft function, but include the non- 
relativistic logarithms only at fixed order \n/5. 

The renormalization group equations and the for¬ 
mula for the resummed cross section are very similar to 
those derived for Drell-Yan production within the SCET 
framework 111 1414116l . with the HH' pair in represent¬ 
ation Ra replacing the colour-singlet Drell-Yan pair or 
electroweak gauge boson. The resummed partonic cross 
section then reads II105II1061 




v-' ; r” JrJE - I OJ I /^\ \e-2rir'; 


(58) 


with 77 = 2flr(Ati,A</)- The summed logarithms are contained in the evolution function 




4M2 

eI ) 


Rh 


X exp 


4(5 (///„///)-5 (/ts,ju/)) 


■ la](p.h,Rs) + 2a‘^’''(yu„7t/) H- la"''''' 


(59) 


and 




f 


wf"(w,yu) 


(60) 


denotes the Laplace-transform of the MS-renormalized 
soft function with respect to s = The 

sum over the final-state representations Ra in the fac¬ 


torization formula (43 1 has disappeared in the colour 
basis 0 , since there is a unique final-state represent¬ 
ation for each term in the sum over i. The resummed 
partonic cross section is then matched to the fixed-order 
cross section at the highest available order (NNLO for 
top quarks, NLO for supersymmetric particles), and in¬ 
tegrated with the parton luminosity. The presence of 
the Coulomb functions introduces some subtleties in the 
convolution with the parton densities, which are dis¬ 
cussed in 111 171 . 


For the definitions of the quantities appearing in (58 1 , 
( |59| ) we refer to 11061 . and mention only briefly the ne¬ 
cessary ingredients for NNLL resummation. The hard 
functions are process-specific and must be obtained 
from the one-loop 2^2 production cross sections dir¬ 
ectly at threshold, separated into irreducible colour (and 
if necessary, spin) representations. The soft function 
is also needed at one-loop, and has been computed for 


arbitrary colour representations in 11051 . The single¬ 
particle and cusp anomalous dimension must be used 
at the two-loop and three-loop order, respectively, and 
are the same as for the Drell-Yan process. The only 
new anomalous dimension is the anomalous dimension 
for the soft function (56 1 . In 11051 the required ex¬ 
pression has been related to the constant coefficient in 
the anomalous dimension of the heavy-heavy current in 
heavy-quark effective theory in the limit where the cusp 
angle goes to infinity. At the two-loop order relevant 
to NNLL resummation, the soft anomalous dimension 
exhibits Casimir scaling, 7 *"^ = Cr^jh.s- From the ex¬ 
pression for the two-loop anomalous dimension of the 
heavy-heavy current 011811119l . one extracts 


r!!.’ = 




2 jt ^ \ 40 

^ + 4 ^ 3 )+ 


( 61 ) 


for the coefficient of (aJiAn))^. This result was con¬ 
firmed by an explicit, independent calculation 01201 . 


4.2. Top-pair invariant mass distribution near 
threshold 

A Coulomb enhancement of the cross section for the 
production of top-quark pairs close to threshold, which 
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has been discussed in Section [3] in the context of 
collisions, can also be observed at hadron colliders, 
since it is possible to produce the top quarks in the 
colour-singlet state. Indeed, at LHC, the cross section 
close to threshold in dominated by the process gg tt 
where the top-quark pair is in the *5o colour-singlet 
state. In contrast to a linear collider, where the physical 
observable is the total cross section as a function of en¬ 
ergy, at the hadron collider one considers the invariant- 
mass distribution of the top-quark pairs. 

The calculation of the cross section within the 
NRQCD framework contains as building blocks the 
hard production cross section for a top-quark pair at 
threshold ifnn wm and the non-relativistic Green 
function governing the dynamics of the would-be to- 
ponium bound-state. In 012111 the NLO formulae 
were derived for quark or gluon initial states and a 
quarkonium in a = 0“^ colour-singlet state, plus 
possibly a parton. The general case, with the heavy 
quark system {QQ) in S-wave singlet/triplet spin state, 
and colour-singlet/octet configuration is given in 01220. 
together with the corresponding results for P-waves. 
The results of Ifnn fmi were presented for stable 
bound states. For wide resonances it is convenient to de¬ 
scribe the bound-state dynamics through a Green func¬ 
tion. 

NLO calculations have been performed in 0111II1121I . 
where slightly different approaches have been applied. 
Whereas in 01 121 the matching has been performed in 
the strict threshold limit where the partonic centre-of- 
mass energy s approaches twice the top-quark mass, the 
complete dependence on s as given in 0121111221 is in¬ 
cluded in 01 111 . Thus, formally, the result of 01121 is 
only valid for top-quark production where the velocity 
of both quarks is small. In the approach of Iml on the 
other hand, the relative velocity still has to be small but 
the combined top-antitop quark system can move with 
high velocity. Furthermore, Ref 01 1 II includes all NLO 
subprocesses, i.e. also those which appear for the first 
time in (9(a^) and performs a soft gluon resummation at 
the NLL order using the Mellin-space approach. 

In Figure[^the invariant-mass distributions for LHC 
with Vs = 14 TeV centre-of-mass energy is shown for 
the sum of all contributing channels and separately for 
the colour-octet and colour-singlet contribution. The 
width of the bands is obtained from varying renormal¬ 
ization and factorization scales in the hard cross section 
as described above. The additional uncertainty from the 
Green function, which is estimated to about 20% for the 
singlet and below 5% for the octet case Hill , is not in¬ 
cluded. 

As expected, for invariant mass M < 2m, the pro- 



M IGeV] 


Figure 10: Invariant-mass distribution for top-quark pair production. 
The width of the bands reflect the scale dependence of the hard scat¬ 
tering parts. Figure from nm . 


duction of tt pairs is dominated by the singlet contri¬ 
bution. However, for M > 2m, one observes a strong 
rise of the octet contributions, in particular of the gluon- 
induced subprocess, which for M > 2m, -H 5 GeV be¬ 
comes even larger than the corresponding singlet contri¬ 
bution. For the colour-octet case the scale dependence 
of the hard scattering amounts to +7%. Considering the 
threshold behaviour as shown in Figure[^it is clear that 
the location of the threshold is entirely governed by the 
behaviour of the colour-singlet (S-wave) contribution. 
Thus, as a matter of principle, determining the location 
of this step experimentally would allow for a top-quark 
mass measurement, which is conceptually very different 
from the one based on the reconstruction of a (coloured) 
single quark in the decay chain t —> Wb. In fact, much 
of the detailed investigations of tt threshold production 
at a linear collider were performed for this situation by 
establishing the relation between the location of the col¬ 
our singlet top-antitop resonance and the top-quark MS- 
mass. 

Considering the threshold region (say up to M,f - 
350 GeV) separately, an integrated cross section of 
15 pb is obtained, which should be compared to 5 pb 
as derived from the NLO predictions using a stable top 
quark and neglecting the binding Coulomb-force cor¬ 
rection. Within this relatively narrow region the en¬ 
hancement amounts to roughly a factor three and a sig¬ 
nificant shift of the threshold. Compared to the total 
cross section for tt production of about 900 pb, the in¬ 
crease is relatively small, about 1%. However, in view 
of the present and future experimental precision these 
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effects should not be ignored. 

4.3. Total top pair production cross section 

The ability of LHC to produce top-quark pairs in 
large numbers has triggered a large amount of theor¬ 
etical work devoted to improving the accuracy of the 
prediction of the total cross section. The first results re¬ 
ferred to approximate NNLO (“NNLOapp”) calculations 
that aimed at including the singular terms in the partonic 
cross section as/? ^ 0 0123141261 . Soft-gluon resumma¬ 
tion of the total cross section with NNLL accuracy was 
completed in 0117111271 with two independent calcula¬ 
tions, one done in the joint soft-Coulomb resummation 
approach described above, the other in the Mellin-space 
soft-gluon resummation formalism. In yet another ap¬ 
proach the total cross section is computed from NNLL- 
resummed or approximated NNLO calculations of cer¬ 
tain differential distributions 11128141301 . 

It should be noted that threshold resummation for the 
top-pair inclusive cross section does not have a clear 
parametric justification. In order to obtain the total had¬ 
ronic cross section, the partonic cross section is con¬ 
voluted with the parton distribution functions. Both at 
Tevatron and LHC, the top-antitop invariant-mass distri¬ 
bution peaks at about 400 GeV, which corresponds (in 
the absence of radiation) to j3 ^ 0.5. The convolution 
of the partonic cross section with the parton luminos¬ 
ity is therefore dominated by the region j3 > 0.3, where 
the threshold approximation is no longer valid. Nev¬ 
ertheless, one often hnds that the threshold expansion 
provides a reasonable approximation even outside its 
domain of validity. At the very least, the approximation 
is better than the one not using this piece of information. 

Much of the ambiguity in resummed calculations has 
been removed through the hxed-order NNLO calcula¬ 
tion II1311I1321 . The present state-of-the-art prediction 
therefore consists of NNLL resummation matched to 
the full NNLO result, and is available in the programs 
Topix^ lll33l (based on 111 171 1 and top-n- 113411 (based 
on 11271 1. Topixs includes soft and Coulomb resumma¬ 
tion and is therefore technically the most complete the¬ 
oretical prediction. In particular, it includes the bound- 
state effects in the threshold region discussed in the pre¬ 
vious subsection. However, the effect of Coulomb re¬ 
summation beyond the terms already included in the 
hxed-order NNLO result is rather small for the total 
cross section. 

In Table we summarize successive approximations 
to the top-quark pair production cross section at Tev¬ 
atron and LHC from NLO to NNLL, where NNLL 


^http://users.ph.tum.de/t31 software/topixs/ 


o-„-[pb] 

Tevatron 

LHC ( ^/s =n TeV) 

LHC ( ^/s =8 TeV) 

NLO 

A AS+0-36-)-0.23 
°-°^-0.75-0.22 

1 1 -bl9.5-b6.8 

1D0.1_2i 2_6.2 

74-27.84-9.2 

z,z,u.^_297_g3 

NNLOapp 

7 i)f:+0.26+0.29 
'•UO_0.34-0.24 

l + 12.3-b7.3 

iUi.i_ii 9_g 7 

Z.OU.U_jg2_gQ 

NNLO 

7 M1-bO.27-1-0.29 
'•Ul_o.37_o.24 

167 6.7-t-7.7 

iU/.l_io7_7 1 

99Q 1+9.24-10.4 

NNLL 

7 1 c-b0.24-)-0.30 
' “^-0.10-0.25 

168 5+ 6.34-7.7 

75_72 

941 a+8.74-10.5 
^^^•^-11.1-9.7 


Table 1: Top-quark pair production cross section at Tevatron and 
LHC for m, = 173.3 GeV, a^iMz) = 0.1171 ± 0.0014, (N)NLO 
MSTW08 PDFs 11351 . The first error represents to theoretical un¬ 
certainty from independent soft/hard/Coulomh scale variations and 
resummation ambiguities in the partonic cross section, the second 
PDF+ff., at 68% CL. 


means resummed matched to full NNLO, generated 
with Topixs 2.0. We note that in production at the Tev¬ 
atron, which is dominated by the quark-antiquark ini¬ 
tial state, the resummation/NNLO effect is signihcant 
(-1-8%, NNLL vs. NLO in the table), and the threshold 
approximation to the full result worked well (NNLO 
vs. NNLOapp). On the other hand, in gluon-gluon ini¬ 
tiated production at the LHC, resummation is a small 
correction (-1-1%) and underestimates the full NNLO 
correction (-1-4%). However, in both cases, resumma¬ 
tion without the full NNLO correction already leads 
to a significant reduction of the theoretical uncertainty 
(8% ^ 3% at Tevatron, 13% ^ 4.5% at LHC, exclud¬ 
ing the PDF -H as error), while still including the full 
NNLO-bNNLL result in the uncertainty estimate. 

Figure 11 summarizes the results for the total top- 
pair cross section at the Tevatron from different theor¬ 
etical calculations. By comparing the black and blue 
bars to the left and right of the vertical dashed line, we 
observe that the predictions from joint soft-Coulomb 
resummation (first bar [black], il33l ') and the Mellin- 
space approach (second bar [blue], 013211 ') are in very 
good agreement, once the full NNLO result is included 
and resolves some of the resummation ambiguities (in 
favour of the former). The NNLOh-NNLL computations 
are in good agreement with the measurement for the ad¬ 
opted top-quark pole mass m, - 173.3 Ge'V, while the 
NLO result would underestimate the measurement sig¬ 
nificantly (compare Table[^. 


The top-quark pair cross section has become a pre¬ 
cisely predicted and measured quantity, which depends 
essentially only on the fundamental parameters m, and 
as, and the parton distributions. Assuming standard 
physics the measurement constrains these parameters. 
The possibility to determine the top-quark mass in a the¬ 
oretically clean way, though less precisely than from re¬ 
construction of top decay products, has been explored. 
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Tevatron, NNLL 

a-[pb] 

8.0 - ! 


6.0 - * I 

5.5 ' 

Figure 11: Top-pair production cross section calculations with NNLL 
accuracy in pp collisions at = 1.96 TeV (Tevatron). To the left of 
the dashed line: NNLOapp + NNLL results. Black: Beneke et al. 
11171 and Topixs 1.0 (m, = 173.3 GeV); Blue: Cacciari et al. 11271 
(m, = 173.3 GeV) Red: Ahrens et al. (IPI, lUt = 173.1 GeV) 
Magenta: Ahrens et al. 11291 (PIM, m, = 173.1 GeV) To the right 
of the dashed line: full NNLO + NNLL results. Black: Update of 
Beneke et al. 11331 and Topixs 2.0 = 173.3 GeV) Blue: Czakon et 

al. 11321 and top++2.0 {m, = 173.3 GeV) Darkgreen: CDF/DO com¬ 
bined measurement quoted from CMS/ATLAS-CONF-2012-149. Er¬ 
ror bars: inner solid - theory uncertainty of the pailonic cross section 
excluding a^, outer dashed - PDF (+ 0 * 5 ) added. PDF set: MSTW2008 
NNLO with OsiMz) = 0.1171 ± 0.0014. 

for instance, in 111 171 1133111361 11371 . The determina¬ 
tion of the strong coupling has been considered in 0138L 
The impact of the top-quark cross section at LHC on 
the gluon distribution in the proton has been investig¬ 
ated first in 11331 . and in more detail in II139I within the 
NNPDF framework. See also II 1401 . 

4.4. Pair production of supersymmetric particles 

The search for the partners of the SM particles pre¬ 
dicted by supersymmetric extensions of the SM is an 
integral part of the LHC physics programme. Present 
exclusions already imply that the masses of the strongly 
interacting squarks and gluinos are most likely in the 
TeV range. For such heavy particles, threshold resum¬ 
mation is expected to be important, and, contrary to the 
case of top quarks, even parametrically relevant, as a lar¬ 
ger fraction is produced close to the partonic threshold 
due to the fall-off of the parton distributions at large mo¬ 
mentum fraction. 

The study of soft-gluon resummation for pair pro¬ 
duction of supersymmetric particles was initiated by 
111 13lll4l1l . The NLL analysis of the squark-antisquark, 
squark-squark, (anti) squark-gluino and gluino-gluino 
final states 01421 finds significant corrections of sev¬ 
eral 10% beyond the fixed-order NLO calculation, es¬ 
pecially for the gluino-gluino final state, which involves 
the largest colour charges. Note, that NLO correc¬ 
tions to gluino bound-state production has been con¬ 


sidered in 0143111441 and gluino pair production close 
to threshold is presented in 11451 . 

The resummation formalism for soft and Coulomb 
gluons was in fact first used to predict squark-antisquark 
production 01061 in the NLL approximation as defined 
in ( [57| ) and then extended to all superparticle pair final 
states in EU, which also generalized the factorization 
formula to a particular case of P-wave production rel¬ 
evant to stop-antistop production. A scenario with su- 
perparticles with substantial decay widths (for instance, 
gluinos decaying further into squarks and quarks) was 
also investigated 01471 . which adds the complication of 
unstable-particle effects and non-resonant production. 

Different from top-pair production, the summation 
of Coulomb effects is important for some final states, 
which becomes apparent in large differences between 
the NLL results of 01421 and 01(^111461 . The size of the 
resummation effects can be quantified by 

jjjatched 

■^NLL - -. \OZ) 

CTNLO 

where ctnll, matched is the resummed result, properly 
matched to the full NLO calculation. The NLL soft- 
and Coulomb-resummed result is shown as solid line in 
Figure [T^ the one without Coulomb summation is the 
dashed NLLj+h line. In squark-antisquark production 
(upper panel), for large squark masses, the resummation 
effect is more than a factor of four larger in the former 
treatment, which highlights the effect of Coulomb at¬ 
traction. No such effect is observed in the squark-squark 
production process (middle panel), due to a cancellation 
between the numerically dominant, repulsive colour- 
sextet channel for same-flavour squark production, and 
an attractive colour-triplet channel in different-flavour 
squark production 01461 . Gluino-gluino production lies 
in between these extreme cases, but shows the largest 
resummation effects as already mentioned above. Fig- 
ure[^also demonstrates that bound-state production (or 
rather, the resonant enhancement in and below the nom¬ 
inal threshold region) provides a significant further en¬ 
hancement of the total cross section, given by the differ¬ 
ence between the dot-dashed and solid lines. 

The main difference between nm and 0106111461 
can be traced to the soft and Coulomb interference 
term at NNLO, which is already included in NLL soft- 
Coulomb resummation, but not in NLL soft-gluon re¬ 
summation alone. The difference should therefore be 
reduced, when soft-gluon without Coulomb resumma¬ 
tion is extended to NNLL as has indeed been confirmed 
by 01481 . NNLL results in both frameworks have mean¬ 
while been presented 0149111500 for all final states, and 
an additional result is available for gluino pairs nni. 
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PP -> qq (s/7 = 14 TeV) 

^NLL 



Figure 12: NLL K-factor for squark-antisquark (top), squark-squark 
(middle), and gluino-gluino (bottom) production at LHC with = 
14 TeV. The plots show ^nll as a function of M = for 

dilferent NLL approximations: NLL (solid blue), NLLnoBS (bound- 
state contributions excluded, dot-dashed purple) and NLLs+h (dashed 
red). See the text for explanation. Figure from QS). 

5. Sommerfeld enhancement of dark matter pair 

annihilation 

In this last section we turn to a potentially import¬ 
ant non-relativistic effect in the pair annihilation of dark 
matter (DM) particles. While the nature and origin of 
dark matter are still unknown it is intriguing that the ob¬ 
served abundance can be explained rather naturally as 
thermal relic of a TeV scale particle with weak inter¬ 
action strength. The relic density is determined by the 
total annihilation cross section, (crv), averaged over the 
velocity distribution of the particles. Since the abund¬ 
ance froze out when the temperature of the Universe 
was about M^/25, where is the DM particle mass, 
the typical velocities v ~ 0.2 are non-relativistic. When 
DM particles annihilate in the present Universe, poten¬ 


tially revealing themselves in cosmic ray signatures, the 
typical velocities are v ~ 10“^, and the annihilation oc¬ 
curs even deeper in the non-relativistic regime. 

For heavy, weakly interacting, TeV-scale dark mat¬ 
ter the exchange of the electroweak gauge (and Higgs) 
bosons between the slowly moving DM particles gen¬ 
erates a Yukawa potential, which leads to an enhance¬ 
ment of ladder diagrams. Unlike the case of the long- 
range Coulomb potential due to gluon or photon ex¬ 
change discussed up to now, the enhancement does not 
increase as aEw/v at very small velocities, but is cut 
off at aEwM^fMw by the mass of the mediator particle. 
Nevertheless, additional particle exchange is not sup¬ 
pressed by the size of the electroweak coupling qew 
when the DM mass is much larger than the mediator 
mass. This so-called Sommerfeld effect can exceed 
the lowest-order annihilation cross section by orders of 
magnitude. The Yukawa potential has only a finite num¬ 
ber of bound states of order OiaEwM^IMw). When the 
dark matter mass is varied, resonant enhancements of 
the DM pair scattering wave-functions occur whenever 
a bound state approaches the threshold, which can yield 
pronounced Sommerfeld enhancements stronger than in 
the Coulomb case. 

The relevance of the Sommerfeld effect was first 
pointed out for (wino- or higgsino-like) neutralino DM 
annihilation into two photons 01521 . and subsequently 
for relic-density calculations 01530 . although it was not 
until 2008, when an anomalous positron excess was 
measured by PAMELA, that Sommerfeld enhanced DM 
models attracted more attention as a mechanism to boost 
the DM annihilation rates 01540 . We note that for heavy 
dark matter with mass far above the electroweak scale, 
the mass splittings between the states of the DM elec¬ 
troweak multiplet are naturally in the GeV or sub-GeV 
range. At freeze-out a multitude of co-annihilation pro¬ 
cesses, including charged states are still active, and 
Sommerfeld enhancements are generic. This occurs 
in particular in the minimal supersymmetric standard 
model (MSSM). Sommerfeld enhancements in cosmic 
ray signatures and in the thermal relic abundance have 
therefore been discussed extensively for relevant MSSM 
scenarios with neutralino DM 01551416l1l . but also for 
generic multi-state dark matter models 0162141641 . 

In the following we summarize the work of 01654 - 
UMl, which aims at improving the calculation of the 
dark-matter annihilation cross section and relic abund¬ 
ance by including the Sommerfeld radiative correc¬ 
tions in the general MSSM, beyond the previously con¬ 
sidered wino- and Higgsino-limit, where the neutralino 
can be an arbitrary admixture of the electroweak wino, 
Higgsino and bino states. As in previous sections a non- 
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Figure 13: Graphical representation of the I I forward-amplitude. 
The square stand for the action of the potential (the sum of ladder dia¬ 
gram) that scatters I into another state I' (J' in the conjugated amp¬ 
litude), which then annihilates into a final state of Standard Model 
particles (and MSSM Higgs bosons, if kinematically allowed). [Fig¬ 
ure: courtesy C. Hellmann] 


relativistic effective theory (of the MSSM) is construc¬ 
ted to separate the short-distance annihilation process 
from the long-distance Sommerfeld effect, which is en¬ 
coded in the matrix elements of local four-fermion oper¬ 
ators. The approach is very similar to the NRQCD treat¬ 
ment of quarkonium annihilation El, except that we 
deal with scattering states of several species of particles 
interacting through the electroweak Yukawa force. 

The presence of several particle species xi entails 
many complications. Since electroweak gauge boson 
exchange may change the two-particle state (for in¬ 
stance, scatter a neutralino pair into a pair of oppos¬ 
itely charged, nearly degenerate, charginos), the short- 
distance annihilation process is described by a matrix 
in the space of two-particle states, which is not diag¬ 
onal. The off-diagonal terms cannot be obtained from 
the tree-level cross sections computed by numerical pro¬ 
grams. The non-relativistic expansion of the short- 
distance annihilation cross section is therefore com¬ 
puted analytically rather than numerically, which is also 
necessary to separate S- and P-wave, and spin-0 and 
spin-1 annihilation contributions, which receive differ¬ 
ent Sommerfeld enhancement factors. The large num¬ 
ber of final and initial states in the MSSM imply that 
several thousand processes have to be computed. The 
long-distance part of the problem then involves the solu¬ 
tion of a matrix-valued Schrodinger equation in the 


space oiXiXj two-particle states. The existence of kin¬ 
ematically closed channels leads to numerical instabil¬ 
ities that requires a new method to determine the scat¬ 
tering wave functions at the origin relevant to the Som¬ 
merfeld enhancement. These issues together with a few 
results are summarized below. For many further relev¬ 
ant details, we refer to 0165111681 . 


5.1. Construction of the annihilation rates 

The process we are interested in is shown schematic¬ 
ally in Figure [T^and is described by the non-relativistic 
MSSM Lagrangian 

^NRMSSM ^ ^ ^ ^ (63) 

Xkin contains the bilinear terms in the two-component 
spinor fields ft and xpj - rjj, that represent the no < 4 
non-relativistic neutralinos (x^) and n+ < 2 charginos 
(xj wdxp, respectively. The entire treatment is restric¬ 
ted to the leading-order terms in the velocity expansion 
in the long-distance part, hence Xkin is simply given by 


Xkin - 

/-I 

n+ 


id, - (mi - otlsp) 


2otlsp 


j=i 


f2 f 


id, - (nij - otlsp) + 


Ini] 


LSP 


(64) 


The non-relativistic energy is measured relative to the 
mass mLSP of the lightest neutralino state. To have a 
consistent power-counting, the mass differences (m,- - 
wiLSp) must formally be considered of order mLSpv^- 
The short-distance annihilation of the chargino and 
neutralino pairs into SM and light Higgs final states is 
reproduced in the effective field theory by local four- 
fermion operators. The leading-order contributions to 
•J-Cann ^Tc given by dimension-6 four-fermion operators, 
that describe leading-order S-wave neutralino and char¬ 
gino scattering processes ^ YeiYej- They read 




d=6 

ann 


V' 1 fXX-^XX (2S + la ) (~XX^XX (2S + \a ) 

2_J Z_J 4 ■^{ei^^2)(c4e3) V MC4C3)(e2<’l) V ’ 

XX^XXS=Q,\ 


(65) 


where are the corresponding short- 

distance coefficients. The explicit form of the 
dimension-6 S-wave operators with spin S = 0,1 is 

3 ( 66 ) 

= xlyXe^XeWXe, ■ (67) 


The first sum in ( |65] l is taken over all neutralino and 
chargino forward-scattering reactions ^ Yc 4 Ye 3 3 

including redundant ones where the particle species at 
the first and second and/or third and fourth position are 
interchanged. Thus in the sector of two-particle states 
with electric charge Q, is a square matrix of 
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Figure 14: Matching of the imaginary part of the forward amplitude. 
[Figure: courtesy C. Hellmann] 


dimension = (No,Ni,N 2 ) — (24,16,4). In order 
to reproduce the short-distance tree-level annihilation 
cross section with Oiy^) accuracy, we add dimension- 
8 operators with two derivatives accounting for P-wave 
and corrections to S-wave annihilation, but also operat¬ 
ors proportional to the mass differences 6m - - 

m(,j)/2, 6m - {m^^ — 171^^12. Of course, the latter have 
non-vanishing coefficients only for the off-diagonal ele¬ 
ments of forward amplitude. 

The short-distance coefficients of all these operators 
are determined by expanding the MSSM amplitudes for 
the process XeiXe 2 X ^ XetXej, with SM and light 
Higgs intermediate states in the relative momenta p, p' 
and mass differences 6m, 6m, and matching them to the 
tree-level matrix element of the four-fermion operators 
for the same incoming and outgoing states to 0(v^), as 
illustrated in Figure[^ For the computation of the neut- 
ralino and chargino inclusive annihilation rates, only the 
absorptive part of the short-distance coefficients are re¬ 
quired according to the optical theorem. At tree-level, 
it is also possible to dehne the rates to exclusive hnal 
states (relevant to indirect detection of DM) in this way. 
We note that the matching is performed for on-shell 
scattering, which implies that p and p' are different for 
the off-diagonal scatterings. Energy conservation im- 

= yfs — M + {6m + 6m ) -H ... , (68) 

= s/s — M — (6m ■+■ 6m ) -H ..., (69) 

where ^fs is the centre-of mass energy of the scattering, 
M - (liti m - (m^j -H meJ/2, m - (m^^ -H 

me,)l2, and p - mmIM. The non-relativistic expansion 
is strictly valid, when sfs - M, 6m and 6m are of order 
M^. 

The actual calculation is complicated by the pres¬ 
ence of many hnal states and interactions in the MSSM. 
Moreover, for reasons discussed in Gsa, the calcula¬ 
tion is most easily done in ‘t Hooft-Feynman gauge, 
which, however, adds a large number of unphysical h- 
nal states containing pseudo-Goldstone Higgs and ghost 


plies that 

t 

2p 

e1 

2p 




Figure 15: Numerical comparison of the non-relativistic approxim¬ 
ation (solid lines) to the tree-level annihilation cross-section times 
relative velocity, crvrei, for (upper figure) and 

X^X 2 H^H~ (lower figure) with the corresponding unexpan¬ 
ded annihilation cross section produced with MadGraph fT69l . The 
second process is dominated by P-wave annihilation. Vrei is given by 
Vj.gj = for tho, XeiXei process. The underlying 

MSSM spectrum is a wino-like neutralino LSP scenario, generated 
with the spectrum calculator SuSpect flTOl . The masses of the 
and^J^ are given by m^o = 2748.92 GeV, m^+ = 2749.13 GeV, and 
m^- = 3073.31 GeV. The mass of the Higgs particles takes the 
value m//± = 167.29 GeV. Figure from 11651 . 


particles. The calculation of all S-wave annihilation 
rates is presented in 01651 . and the coefficients of the 
0{v^) terms in 01661 . 

The non-relativistic expansion to 0{v^) is usually a 
very good approximation. For the diagonal annihila¬ 
tion rates, the result from the above calculation can be 
compared to the unexpanded tree-level cross section. A 
comparison for two cases is shown in Figure [TS] Good 
agreement is found up to relative velocities of about 
0.5, which is sufficient for all practical purposes. The 
analytic expressions from insmsa cover all relev¬ 
ant processes, when only neutralino and chargino (co-) 
annihilation processes are relevant. MSSM scenarios 
with sfermion-neutralino co-annihilation are presently 
not covered, neither are resonant annihilation models, 
since in this case the annihilation process is no-longer 
short-distance, and cannot be expanded in the velocity. 
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5.2. Computation of the Sommetfeld enhancement 

Once the short-distance coefficients are at hand, the 
Sommerfeld-enhancement is computed from the matrix 


elements of four-fermion operators such as 

The first step consists in deriving the potentials in the 

expression 


Aot = - Xi \ (70) 

XX^XX ^ 

from all XexXei Xe^Xe, neutral, single-charged and double-charged scattering reactions through electroweak gauge 
boson, photon and Higgs boson exchange. Assigning mass m^ to the exchanged particle, the potential in coordinate 
space has the form 

- (^eie2e4e3^<t4Q'l^a3Q'2 


(71) 


which contains spin-independent and spin-dependent 
terms. The indices a, are contracted with the (unwrit¬ 
ten) spin indices of the field operators in ( [70| ), and the 
total spin operator ^ is ^a 4 ai,a,a 2 = ^aiaj' 26 a^a 2 + 
Since the total spin is not changed at 
this order, the two-particle states XeiXe 2 Xe^Xe^ un¬ 
dergoing potential interactions can be decomposed into 
partial-wave states with defined spin S = 0,1, 
and the Sommerfeld factors can be determined separ¬ 
ately for each spin. The MSSM neutralino-chargino po¬ 


tentials are given in 0167L 

The resummation of ladder diagrams carrying the 
Sommerfeld enhancement is achieved by including the 
interaction ( |70| ) into the unperturbed Lagrangian. The 
task is then to evaluate the matrix elements of the anni¬ 
hilation operators in this theory. For the general treat¬ 
ment we refer to II 16711 and consider here the example of 
the '5o operator (jh^. Its matrix element can be para¬ 
metrized in the form 


(’Ci^j\^{e,7Me2e,}C^0)\XiXj) = <XiX jlXe.Xe^l^} {^IxllXe, \XiX j) 




(72) 


where (in general) is the -component of the 

scattering wave function for an incomingstate with 
centre-of-mass energy yfs, orbital quantum number L 
and total spin S , evaluated for zero relative distance and 
normalized to the free scattering solution. In the ab¬ 
sence of potential interactions, the tree-level matrix ele¬ 
ment of the four-fermion operators is obtained by repla¬ 
cing —> dej 6 e^j. It is convenient to introduce 

a notation where a single index refers to a two-particle 
state (for instance, e = e\e2,i - if) rather than two 
indices for the individual particle species that make up 
this state. The Sommerfeld factor is defined as the anni¬ 
hilation cross section including the potential interaction 
relative to the tree cross section. For a given partial- 


wave contribution to the annihilation cross section, it is 
defined as 






(73) 


Here fee’ is related to the absorptive part of fe^^^^- A 
sum over all two-particle states e, e' is implicit in (j7^. 

The annihilation cross section including the Sommer¬ 
feld corrections is constructed by multiplying each spin 
and partial-wave component by its respective Sommer¬ 
feld factor fT^ . The tree-level cross section expanded to 
(9(v^) is recovered by setting 5,■[/(^‘^^'Ly)] ^ 1. Defin¬ 
ing p? = 2pi (yfs- Mi) -H 0 {pf ), the relative momentum 
squared of the two annihilating neutralinos, and M, and 
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fii, the total and reduced mass, respectively, of the two- takes the form 
particle system, the master formula for the cross section 


tr-- • v,ei = 5,[/,(‘5o)] ffr^^CSo) + sm^s i)] 3 i) 

i2 , 


light 

^ (si[g.eso)] ^r"’^cso)+sA8.es,)] + 5 , 

\fepj) 


fCPi) 

m2 




+ s. 




ffn’^i^pj)). 


(74) 


The precise definition of the annihilation matrices ap¬ 
pearing in this equation is given in 11671 . 

It is well-known that the matrix elements of four- 
fermion operators can be expressed (at this order in the 
non-relativistic expansion) in terms of wave functions 
at the origin. In the present case - {‘PsiQ)]],, with 
if/Eir) the matrix-valued solution of the Schrodinger 
equation 


V2 

»Jlsp 


.ab 


V‘‘\r) 


[^E{r)]hi = 0 . 


(75) 


Here E - yfs — IrriESp and 

V“>’(r) = y“*(r) + 5“'’ [m, - 2mLSp] ■ (76) 

We also use the velocity variable v defined by 

E = otlspv^ ■ (77) 


The solution method proposed in 117111 (and reviewed 
in 11671 in the present notation) relates the Sommerfeld 
enhancement to the asymptotic behaviour of the radial 
partial-wave function u{r) for certain boundary condi¬ 
tions at r = 0. For the S-wave case, one finds 


[<^£(0)]:,- = [U-\r^)]ei, (78) 

where U is the matrix 

Uab(roo) = {[u(roo)]ab - ika[u(roo)]ab), (79) 

and ^2 = nii^spiE + ie - [M^, - 2mLSp]) for a two-particle 
state a with mass M^. The Schrodinger equation must 
be solved from r = 0 to a value r^x, large enough such 
that Uabifoc) is close enough to its asymptotic value. 

The procedure described here works well, when all N 
states included in the multi-channel Schrodinger equa¬ 
tion are degenerate to a high degree, but fails other¬ 
wise. This is illustrated in Figure which shows 



Figure 16: S-wave spin-0 Sommerfeld factor for v = 0.012 in a wino- 
like model as described in the text. The light-grey (red) dashed curve 
shows S{Xoo), when only the x^[X^[ two-particle states in 
the Schrodinger equation and the annihilation process are kept and 
the asymptotic regime is reached for .Voo > 50. The dark-grey (blue) 
solid curve shows the result when the x^iX 2 included. In this 

case, the evaluation fails for Xco > 2 and no reliable result is obtained. 
Figure from GH). 


the Sommerfeld factor for the x'lX, *5o annihilation 
cross section as function of x™ = nruspi'^oo for a 
MSSM parameter point where the lightest neutralino 
(LSP) is wino-like. The relevant masses are otlsp = 
MyO — 2749.4 GeV, OTy+ = 2749.61 GeV and ntyO — 
2950.25 GeV. The light-grey (red), dashed curve shows 
S (Xoo) when only the two highly degenerate 
x\Xi two-particle states are included in the Schrodinger 
equation and the annihilation process. The velocity is 
chosen v = 0.012, slightly below the threshold for the 
XiXi state. After a rapid variation with a peak structure, 
the Sommerfeld factor reaches a plateau and for Xoo > 
50 stays at the constant value S ( 00 ) ss 199.59. When the 
X^X 2 added to the problem, the Schrodinger sys¬ 

tem is extended to a 3 x 3 matrix. Since the new state is 
200 GeV heavier and moreover rather weakly coupled 
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to the two lowest, nearly degenerate wino states, it 
should have little effect on the value of the Sotntnerfeld 
factor. However, now the numerical solution fails when 
Xoo is slightly larger than 1, as seen from the dark-grey 
(blue) solid curve in Figure which drops to 0 after a 
few spikes. It is not possible to reach the plateau, where 
S (Xoo) stabilizes. 

The numerical instability originates from the pres¬ 
ence of kinematically closed two-particle state chan¬ 
nels, here the x^X 2 state. The solution [u{x)]t,i for the 
closed channel involves an exponentially growing com¬ 
ponent proportional to e**'' where = inLSpiMb - 
[2mLSP + rwLSpt'^]). The off-diagonal potentials V^^ir) 
couple the different channels and the open-channel solu¬ 
tions [M;(x)]a, inherit the exponential growth from the 
closed channels. For the two-LSP channel, expo¬ 
nential growth occurs when at least one of the included 
kinematically closed channels b satisfies 

Mb - [2mLSp + 'WLSpv^] > —^ ■ (80) 

niLSP 

Since typically otlsp » Mew for the dark-matter scen¬ 
arios of interest, this condition is easily satisfied unless 
all two-particle states included in the computation are 
degenerate within a few GeV or less. In consequence 
the formally linearly independent solutions [u]ai degen¬ 
erate and the matrix Uai becomes ill-conditioned with 
exponentially growing entries in the rows correspond¬ 
ing to open channels a. The matrix inversion required 
for ( [78] l can no longer be done in practice for Too large 
enough such that the asymptotic regime is reached, 
which causes the instability seen in Figure [T6| 

The solution to the problem provided in 01671 
is based on an adaptation of the modification of 
the variable-phase method for the Schrodinger prob¬ 
lem 11721 . The idea is to write the solution as a linear 
combination of the linearly independent Bessel function 
solutions fa{x), ga(x) of the free Schrodinger problem: 

[u{x)]ai = fa(x)aai(x) - ga(x)/Sai(x) (8 1 ) 

(no sum over a). Defining Uai = aailga, it can be shown 
that the matrix U defined in ( [7^ is asymptotically re¬ 
lated to Oai by 

U,i(x) e''^“"d«,(x). (82) 

One then derives a differential equation system for the 
matrix q'“^(x) from which follows without hav¬ 

ing to explicitly invert the matrix U. 

Leaving aside limitations related to the CPU time 
needed to solve a system of many coupled differen¬ 
tial equations, this method allows to compute the Som- 
merfeld factors reliably also when many non-degenerate 


two-particle channels are present. This makes the Som- 
merfeld enhancement accessible in a larger part of 
the MSSM parameter space, away from the wino or 
Higgsino limit, where the Sommerfeld effect is less dra¬ 
matic but potentially still a large radiative correction. 
To speed up the numerical calculation. Ref. 11671 also 
discusses an approximation to the treatment of heavier 
channels, and a relation between the 0 {v^) suppressed 
S-wave operator matrix elements and the leading-power 
ones. 

5.3. Results 

We briefly discuss a selection of results from II 1681 to 
which we refer for further details. 

5.3.1. Wino-like lightest neutmlino 

Wino-like x^ dark matter arranges into an approxim¬ 
ate SU(2)/, fermion triplet together with the two char- 
gino states A phenomenological MSSM (pMSSM) 
scenario with wino-like x^l is provided by the SUSY 
spectrum with model ID 2392587 in III73I . A meas¬ 
ure for the wino fraction of a given neutralino LSP state 
is the square of the modulus of the neutralino mixing- 
matrix entry Zbf 2 i- For pMSSM scenario 2392587 the 
X^ constitutes a rather pure wino, |Za^ 2 iP = 0.999, with 
a mass niESp = - 1650.664 GeV. The mass of the 

chargino partner is given by - 1650.819 GeV, 
such that 6 m - m^+ - m^o - 0.155 GeV. For comparison 
purposes these values are taken without any modifica¬ 
tion from the spectrum card provided by II173L where 
the mass parameters refer to the DR-scheme. Eventu¬ 
ally, the analysis should be done with one-loop renor¬ 
malized on-shell masses, since the Sommerfeld effect is 
sensitive to the precise value of the mass splitting. 

In Figure we plot the ratio (cr®^v)/(crP‘^'^v) of an¬ 
nihilation rates including long-range interactions, cr^^v, 
over the perturbative tree-level result, for the 

two-particle states x^iX^ ™d XiXi in neutral sector 
of the model as a function of the velocity vlsp of the 
incoming ;^^j’s in their centre-of-mass frame. The en¬ 
hancement peaks in the vicinity of the threshold of the 
heavier neutral state XiX'i t'LSP - 0.014. Well below 
this threshold, the enhancement for the XiX'i system is 
velocity-independent and of (9(10). 

The quantity that enters the Boltzmann equation for 
the neutralino number density is the thermally aver¬ 
aged effective annihilation rate summed over all co¬ 
annihilation channels, {(Teffv). The upper panel of Fig¬ 
ure [T^ shows (cTeffv) as a function of the inverse scaled 
temperature x = m^o/T. As the temperature decreases, 
the Sommerfeld enhancement increases and reaches two 
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Figure 17: The enhancement of the;^j;^j ?ir\.dx\x\ annihilation cross 
sections for Snowmass model 2392587 relative to the perturbative 
tree-level rate, (cr^^v)/(crP®’'*v). The solid lines refer to the calcula¬ 
tion of the Sommerfeld-enhanced rates with off-diagonal entries in the 
annihilation matrices included. The dashed curves show the enhance¬ 
ment with respect to the perturbative cross sections when off-diagonal 
annihilation rates are not considered. The dotted curve labelled “pure- 
Coulomb enh.” shows the enhancement from photon exchange only 
in ihQx\x\ channel. Figure from 11681 . 


orders of magnitude. Around x > 10"^ the number dens¬ 
ities of the xf are so strongly Boltzmann-suppressed 
with respect to the number density despite the small 
mass splitting that the rates of the charginos basically 
play no role in the effective rate (cTeff v), which is then es¬ 
sentially given by annihilations. After decoup¬ 
ling, (crincluding the Sommerfeld enhancements 
becomes constant, which we can infer from the constant 
enhancement factor for thex^X system for very low ve¬ 
locities shown in Figure The difference between the 
dashed and solid lines in Figure and the upper part 
of Figure[T^demonstrates that it is important to include 
correctly the off-diagonal annihilation reactions. 

Since the lightest neutralino is almost a pure wino, 
the qualitative features of model 2392587 are similar 
to the pure-wino studied already in II153L There are, 
however, quantitative differences. A pure-wino model 
contains only an SU(2)z, triplet of x states in addition 
to the SM particles, while the MSSM model 2392587 
features non-decoupled sfermion states at the 2-3 TeV 
scale with non-vanishing couplings of the Xi and xf to 
sfermions and to the (heavier) Higgs states, which re¬ 
duce some of the annihilation rates relative to the pure- 
wino dark matter case. The difference of the thermally 
averaged effective annihilation rate for the two models, 
chosen to have the same x^ mass, is shown in the lower 
panel of Figure [T^ 

The 0{1Q^) effect seen in Figure [T^ at large x is not 
relevant to the dark-matter relic density computation, 
since freeze-out occurs already at x ~ 20. Nevertheless, 



X = m^jT 



Figure 18: Upper panel: The thermally averaged effective anni¬ 
hilation rate (cTeffv) as a function of the scaled inverse temperature 
X = m^lT for Snowmass model 2392587. The two upper (red) curves 
correspond to the Sommerfeld-enhanced annihilation cross sections 
including (solid line) or neglecting (dashed line) the off-diagonals in 
the annihilation matrices. The lower (blue) curve represents (<Teffv) 
obtained from perturbative (tree-level) cross sections. Lower panel: 
Comparison of (<Teff v) for model 2392587 (solid) with the correspond¬ 
ing pure-wino scenario (dot-dashed). Figures from EU. 


the abundance is significantly modified. It is custom¬ 
ary to solve the Boltzmann equation that determines the 
relic density for the yield Y - nj s, defined as the ratio 
of the number density n of all co-annihilating particle 
species divided by the entropy density s in the cosmic 
co-moving frame. Figure [T^shows the ratio of the yield 
Y calculated from Sommerfeld-enhanced cross sections 
in both the pMSSM and the pure-wino model to the 
corresponding results using perturbative cross sections, 
Tpert, as a function of x. Around x ~ 20 the yields in¬ 
cluding Sommerfeld enhancements start to depart from 
the corresponding perturbative results. The enhanced 
annihilation rates delay the freeze-out of interactions, 
which leads to a reduction of the yield Y compared to 
the perturbative result Tpert- The most drastic reduction 
in T/Tpert occurs between x ~ 20 and x ~ 10^. In this 
region the enhancement factors on the cross sections are 
of (9(10) (and not yet (9(10^) as for very large x). Even¬ 
tually, for X > 10^ the particle abundances in both the 
perturbative and Sommerfeld-enhanced calculation are 
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Figure 19: The ratios of the yield Yj Fpert as a function ofx = m o/T, 

X] 

where Y is calculated including the Sommerfeld enhancement on the 
XX annihilation rates while Fpert just uses the perturbative ones. The 
solid (blue) and dashed (black) curves give the results for the Snow- 
mass model 2392587 including and neglecting off-diagonal annihil¬ 
ation rates, respectively. The dot-dashed (red) curve corresponds to 
F/ Fpert(v) in the pure-wino model. Figure from 11681 . 


frozen in, and the fraction F/Fpert remains constant. In 
case of the wino-like model we hnd that the relic densit¬ 
ies calculated from the yield today read QP^'^h^ =0.112 
and - 0.066. Hence, the Sommerfeld effect 

leads to a reduction of the calculated relic abundance of 
around 40% in this model. Neglecting the off-diagonal 
annihilations in the calculation would underestimate the 
effect considerably (dashed curve in Figure [T9|. 

5.3.2. Higgsino-to-wino trajectory 

The formalism described above allows us to con¬ 
sider lightest neutralino states, which are arbitrary ad¬ 
mixtures of the electroweak gauginos and Higgsinos. 
Ref. B168I defines a “trajectory” in the and M 2 MSSM 
parameter space, which interpolates between an almost 
pure-Higgsino (M 2 » fi) and almost pure wino (/r » 
M 2 ). A large bino fraction is excluded by choosing 
Ml = IOM 2 . The trajectory is then chosen such that the 
perturbative relic density computed with the program 
DarkSUSY |174| agrees with the most accurate determ¬ 
ination obtained from the combination of PLANCK, 
WMAP, BAO and high resolution CMB data, = 

0.1187 + 0.0017. The position of 13 models on this tra¬ 
jectory in the p - M 2 plane is shown in Figure]^ For 
details on the other MSSM parameters and the one-loop 
mass renormalization scheme, see 11681 . The models 
can be categorized as higgsino-like with a wino fraction 
of;Fi below 10% but a higgsino fraction \Zyj 2 \ pH-|Ziv4i|^ 
above 0.9 (models 1 - 6), mixed wino-higgsino 
where both the wino and the higgsino fraction lie within 



Figure 20: The fi-M 2 plane with the 13 models defining the higgsino- 
to-wino trajectory, indicated with diamonds. All trajectory models lie 
on the iso-contour for constant relic density = 0.1187 calcu¬ 

lated with DarkSUSY. As reference we also show the iso-contours of 
constant relic densities = 0.095 (lowermost contour-line) 0.15 

and 0.3 (uppermost iso-contour). Figure from (EU. 


0.1 - 0.9 (models 7 - 9), or predominantly wino-like 
with wino fraction above 0.9 (models 10 - 13). 

For each of these models the relic density is com¬ 
puted perturbatively as well as with the Sommerfeld 
enhancement included in the non-relativistic effective 
theory approach discussed above. The result of this 
study is shown in Figure For the higgsino-like and 
mixed models, the perturbative relic densities 
agree very well with the ones calculated with Dark¬ 
SUSY for the same set of input parameters. For the 
wino-like models a difference of up to 8% is observed. 
The Sommerfeld effect is measured by the height of 
the solid-hatched (red) bars, which give relat¬ 
ive to the full height, which represents The 

difference between the two is generally below 10% for 
the most Higgsino-like model 1 - 5, in agreement with 
the hndings for pure-Higgsino models discussed in the 
context of Minimal Dark Matter 11621 . The effect be¬ 
comes more and more significant as the wino fraction 
increases, reaches a maximum for model 10, and then 
decreases again. For model 10, - 0.394, 

which implies that the relic density is overestimated 
by a factor 2.5, when the Sommerfeld effect is neg¬ 
lected. The occurrence of the maximal effect at para¬ 
meter point 10, which features a lightest neutralino mass 
'«LSP = 2320.986 GeV, can be attributed to the above- 
mentioned zero-energy resonance 11521 in the Yukawa 
potential of the;FiYi annihilation channel. 

These results demonstrate that it will be necessary 
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Figure 21: Relic densities Hh- for models 1 - 13 on the higgsino- 
to-wino trajectory. The charts with dotted (black) hatching are the 
perturbative results Bars with dashed (blue) and solid-line 

(red) hatching refer to a calculation with Sommerfeld-enhanced cross 
sections neglecting and properly including olf-diagonal rates, respect¬ 
ively. The grey shaded band comprises Clh^ values within 5% around 
the mean experimental value = 0.1187. The latter value is 

indicated by the black horizontal line and agrees with the DarkSUSY 
result for all 13 MSSM models on the trajectory. Figure from EH. 


to systematically include the Sommerfeld effect when 
MSSM parameter space constraints on heavy-neutralino 
dark matter from direct and indirect searches as well 
as from collider physics are combined with the require¬ 
ment to reproduce, or at least not exceed, the observed 
abundance of dark matter. The formalism and tools de¬ 
veloped in 0165H168ll make it possible to investigate the 
parameter space of the general MSSM, and to identify 
regions where the Sommerfeld effect is not necessarily 
as pronounced as in the previously studied wino limit 
but still constitutes the dominant radiative correction. 


6. Summary 

Non-relativistic physics and its effective description 
play an important role in many different areas of mod¬ 
ern particle physics both in the electroweak and in the 
strong sector. The field of application reaches from pre¬ 
cision determination of SM parameters, to the calcula¬ 
tion of cross sections near thresholds for lepton or had¬ 
ron collider reactions to high order in perturbation the¬ 
ory, and to the description of DM annihilation processes 
in the early Universe. 

In this article we reviewed the production of heavy 
particles close to threshold within the framework of per¬ 
turbative QCD. Special emphasis was put on top-quark 
pair production both at a future electron-positron col¬ 
lider but also at hadron colliders like the Tevatron and 
the LHC. We explained in detail the required methods 


and described the construction of the effective Lagrangi- 
ans for the effective theories NRQCD and PNRQCD. 

The envisaged experimental accuracy for the meas¬ 
urement of the total cross section cr(e+e“ —> it) requires 
perturbative calculations up to the third order. We dis¬ 
cussed several ingredients in detail and described the 
construction of the cross section. Due to the broad spec¬ 
trum of required techniques, which range from three- 
loop vertex corrections in full QCD to the application 
of non-relativistic perturbation theory to third order in¬ 
cluding the corresponding ultrasoft effects, this project, 
now completed, can be considered as a benchmark cal¬ 
culation in the area of perturbative quantum field theory. 

The hadroproduction of top quarks discussed in this 
review relies to a large extent on a factorization formula 
which separates hard and soft scales and the contribu¬ 
tions from the potential region. It is utilized to perform a 
simultaneous resummation of Coulomb effects and soft- 
gluon radiation, which, in combination with the fixed- 
order NNLO result, leads to precise predictions for the 
total cross section. The same formalism is also applied 
to the pair production of supersymmetric particles. 

Pair annihilation of heavy, weakly interacting dark 
matter particles opens a new and fascinating area of 
non-relativistic physics. Non-relativistic enhancements 
of the annihilation rate can be very large despite the 
fact that the force is generated by the exchange of elec¬ 
troweak gauge bosons. With TeV scale dark matter 
particles, degeneracies within the electroweak multiplet 
of DM are generic, leading to a complicated multi¬ 
channel Schrodinger problem. In this article we de¬ 
scribed the formulation and solution of this problem 
in the non-relativistic effective field theory framework 
with applications to neutralino dark matter of the min¬ 
imal supersymmetric standard model. 
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